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MONITORING MOLECULAR INTERACTIONS USING 
PHOTON ARRIVAL-TIME INTERVAL DISTRIBUTION ANALYSIS 

STATEMENT REGARDING FEDERAL FUNDING 
[0001] Tbis invention was made with U.S. Govenment support under U.S. 
Department of Energy Contract No. DE-AC03-76SF00098. The U.S. Government 
has certain rights in this invention. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0002] The present invention relates generally to fluorescence correlation 
spectroscopy (FCS). Particularly, the present invention relates to improved data 
analysis of photon arrival and photon count data typically supplied to FCS analysis 
by utilizing Photon Arrival time Interval Distribution (PAID) analysis techniques. 

2. Description of Related Art 

[0003] The publications and other reference materials referred to herem to describe 
the background of the invention and to provide additional detail regarding its 
practice are hereby incorporated by reference. For convenience, the reference 
materials are numerically referenced and grouped in the appended bibhography. 
The contents of these publications and other reference materials are hereby 
incorporated by reference. 

[0004] With the availability of the complete sequences of genomes of several 
organisms, it is critical to determine the biological function of the proteins coded by 
those genomes. Analysis of protein-protein interactions is important for this process 
smce it can produce protein-protein interaction maps that place each protein in its 
cellular context, fix)m which it is hoped to infer the protein's function. [1, 2]. 
Several existing methods that monitor protein-protein mteractions are: conventional 
and modified yeast two-hybrid systems along with reconstitution systems, phage 
display, fluorescence resonance energy transfer (FRET) methods, mass 
spectrometry, protein chips, and evanescent wave methods [1-5]. 
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[0005] Fluorescence correlation spectroscopy (FCS) and related single-molecule 
methods are important tools for the in vitro analysis of maoromolecular interactions, 
and are potentially useful for in vivo analysis [6. 7]. FCS-related methods can deteci 
these interactions in a distance-independent feshion. unlike FRET. [8]. FCS-related 
methods detect macromolecular interactions by monitoring fluor^cence fluctuations 
that result when fluorescent molecules difiuse or flow across a tightly focused laser 
excitation volume (femtoliter confocal detection volume). At concentrations less 
than 1 nM, the average molecular occupancy of the detection volmne is smaUer tiran 
one. allowing the detection of photon-bursts genemted by single molecules. This is 
commonly referred to as the "low occupancy" regime. A "photon burst" is the set of 
all photons detected from a single molecule during its transit through the coufocal 
detection volume. Analysis of these photon-bursts has been used to measure the 
distribution of molecular properties, such as fluorescence lifetime, polarization 
anisotropy, and fluorescence resonance energy tiansfer (FRET) [9-12]. At 
concentrations between 1 nM-100 nM, the molecular occupancy is stiU low enough 
to be sensitive to the addition or subtraction of one molecule within that volume. 
This is referred to as the "intermediate occupancy" regime. Although it is not 
possible to separate the photons into bursts from single molecules, the resulting 
fluctuating fluorescence signal contains dynamic information about several 
processes such as tianslational diffusion [13], rotational difl&ision [14], intersystem 
crossing to triplet slates [15], and photobleaching [16]. At concentrations grater 
than 100 nM, many molecules occupy the detection volume and the fluctuations are 
averaged out. This is known as the "high occupancy" regime. The primary 
drawback of using FCS-related mediods for monitoring macromolecular interactions 
is the dynamic range over which binding can be detected. These methods are most 
sensitive in the nM concentration regime, whereas binding constants of protein- 
protein interactions often correspond to higher concentrations (this limitation is 
partially offeet by the abihty of FCS-related methods to detect small 
subpopulations). 

10006] Fluorescence bursts or fluctuations are ideally suited to the study of 
macromolecular interactions. Fig. 1 shows how the properties of the sample 
translate into features of the fluorescjence signal for a detection volume with low 
occupancy. Macromolecular interactions such as homo-dimerization and 
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aggregation can be measured using single-channel methods. Fig. lA depicts a 
single-channel measurement on a sample containing a mixture of monomere 
canying one yellow fluorescent label and tetrameis canying four yellow fluorescent 
labels in sohition. The laser excitation profile (shown in green) and the detection 
pinhole define the effective detection volume. As these molecules diffuse in and out 
of the laser excitation profile, bnrsts of fluorescence photons are detected, shown as 
an intensity time trace to the right The three basic characteristics of a single- 
channel photon burst are: (1) the brightness of the bursts {blue arrov^s), which is 
proportional to the number of fluorescence labels detected, (2) the duration of the 
burst (fed arrows), which is related to die difRision time of the molecule across the 
laser beam, and (3) the time between bursts of the same species (green arrows) 
which is inversely proportional to the concentration of that species. The same 
characteristics apply for the fluctuation analysis used at higher concentrations 
(intermediate occupancy). In general, the fluorescence signal from an interacting 
pair of molecules or an' aggregate has different characteristics than that from a free 
single molecule. The complex or aggregate has a larger hydrodynamic radius, 
which results in a longer diffusion time. It also has more labels than fi«e molecules, 
wliich results in increased brightness of the bursts (fluorescence quenching and 
incomplete labeling are ignored at this stage). To most effectively detect binding or 
aggregation, a data analysis scheme that measures all these properties at the same 
needs to be developed. 

10007] For interactions between macromolecules of different types (for example 
hetero-dimerization of two proteins), extending the analysis to two channels 
improves the sensitivity over one channel analysis [8]. The molecules of one type 
are labeled with one color (eg. yellow), and the molecules of the other type are 
labeled with another color (eg. red.) A complex of the two types of molecules has 
bofli labels. This is the situation shown to flie left of Fig. IB. Signal fixan the two 
fluorophores is separated spectrally onto two detector channels, yeUow and red. In 
addition to the ways described for the single-channel case, the binding of two 
molecules labeled with the yellow and red fluorophores can be indicated by the 
detection of simultaneous photon bursts on both channels (ormge arrow in Fig. IB). 
This coincident detection indicates that both fluorescence labels are present, and 
thus the two molecules are bound. 
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[0008] Slnglenrhaimel data reduction and analysis: The task of the data analysis 
perfonned on these fluorescence signals is to extract these parameters using all of 
the infonnation possible, and to present an interpretable graphical representation of 
the data that summarizes the relevant infonnation in the data set. Several methods 
summarized below have been developed to perform this task. All of these methods 
are able to handle vast amounts of data by "reducing" the data to a one- or multi- 
dimensional histogram, from which the information on the diffiasing species is 
extracted. The trick is to reduce the data, but not too much. As much information as 
practically possible should be retained to characterize the sample. 
[0009J Fluorescence Correlation Spectroscopy (FCS): FCS analyzes 
fluorescence fluctuations through the use of the correlation function [13 17 18] 
Correlation functions calculated from the intensity signal reveal the time scall and 
amplitude of various molecular processes, but do not reveal the biightaess of each 
source. In single-chamiel applications, macmmolecular interactions can be detected 
by monitoring the change in diffusion time resulting fiom the interaction of two 
molecules. However, binding often does not produce a large change in diffirsion 
time: for a sphere, doubling the hydrodynamic volume (for instance by binding two 
equally sized subunits) produces only a 26% increase in diffusion time (since the 
diffusion time scales with the hydrodynamic radius, which roughly scales 
as[molecular weight]^'). Therefore, a large change in size is required to measure 
binding using diflusion constants. A further complication is that the shape of the 
bound molecules is also important; there is no general relationship between diffusion 
tune and the number of subunits. For example, a short, rod-like dsDNA fiagment 
will diffuse more slowly than a globular protein with the same volume (compare the 
diffusion constant calculations for a sphere with those for a rod in [19]). 
[OOlOJ Brightness is a reporter of binding events (ignored by FCS) that can in fact 
be more sensitive than the diffusion time. If two interacting macromolecules are 
both labeled, the brightixess of the interacting complex is double the brightiiess of 
the individual subunits. provided that the quantum yield does not change (which is 
not always the case.) Brightaess has the advantage that the shape of the molecule 
does not affect it, unlike the diffusion time. Several methods have been developed 
to use this infonnation. 
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[0011] Moment Analysis of Fluorescence Intensity Distribution (MAFID) and 
higher order correlation amplitudes: Moments of the photon counting histogram 
can be used to monitor occupancy and brightness of labeled macromolecules [20, 
21]. By comparing the values of the mean (first order moment), the variance (second 
order moment), and the third order moment, values for fho occupancy and brightness 
of two species can be extmcted. In this way, macromolecular interactions can be 
monitored by taking advantage of the change in molecular brightness when labeled 
molecules interact. Another method discussed in [22] uses the amplitudes of higher 
order correlations to extract the occupancy and brightness, but turns out to be 
equivalent. 

[0012] Photon Counting Histogram (PCH) and Fluorescence Intensity 
Distribntion Analysis (FIDA): Rather than calculating the moments of the photon 
counting histogram as described above, it is possible to fit the histogram directly, 
thereby usin§ more information to extract brightness and occupancy [23, 24], In this 
way, sub-populations with different brightness can be sq>arated [25]. The PCH and 
FIDA methods differ mainly in their treatment of the shape of the detection volmne. 
PCH has been used to monitor ligand-protein binding equlibria [26], to probe the 
stoichiometry of protein complexes [27], and to study ohgonucleotide-polymer 
interactions [28]. FIDA has been used to probe receptor-ligand interactions in a 
format compa,tible with ultra-high throughput screening [29, 30]. 
[0013] The PCH and FIDA methods contain information about the brightness and 
occupancy of fluorescent species, but lack the information on dynamics contained in 
the correlation function. For a sample with a single species, it is possible to perform 
PCS and PCH or FIDA on the same data set to extract both tbe brightness and 
diffusion time [31], However, if there are multiple species, each with a different 
diffusion time and brightness, there is no direct way to relate each difiRision time 
found to its corresponding brightness. A method that simultaneously tracks 
diffusion time and brightness is necessary to address such heterogeneous samples. 
[0014] Fluorescence Intensity Multiple Distribution Analysis (FEVEDA): By 
using a series of photon count histograms with multiple time bin widths, it is 
possible to obtain the same temporal information as FCS while gaining the 
information on brightness [32, 79]. This is because the shape of the photon count 
histogram is affected by the fluctuations that occur on the time scale of the time bin 
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Width. TMs method is termed Fluorescence Intensity Multiple Distribution Analysis 
(FIMDA). Macromolecular interactions can be tracked using FIMDA by monitoring 
brightness and difiusion time simultaneously. 

[0015J The type of information that is available from each of the previously 
described methods is as follows: FCS extracts concentration and diffiision time (and 
other temporal dynamics), but not brightness; MAFID. PCH. and FIDA extract 
concentration and brightness, but not diffusion time; and FIMDA extracts 
concentration, brightness, and diffusion time (and other temporal dynamics). 
I0016J MuItiplcscbanneJ data reduction and analysis: Coincident detection of 
two fluoiophores of different colors is a more sensitive indicator of binding events 
than brightness or difiusion time used in single^hannel studies [8]. This is because: 
(1) coincident bursts are only detected when two molecules are associated, (2) it is 
less sensitive to quenching of fluorescence, and (3) coincident detection in two 
channels benefits from the properties of ratiometric measmement If two mteracting 
molecules are labeled with the same .fluorophore. it is necessary to detect distinct 
subpopulations with a factor of 2 difierence in brightness. If they are instead labeled 
with different-color fluorophores, the experiment is reduced to a simple "yes or no" 
question. A signal in each channel indicates the presence of the corresponding 
species. A simultaneous signal in both channels indicates a complex (1:1 ratio 
between channels), and a signal on only one channel indicates a free molecule (1:0 
or 0:1 ratio between channels; random coincidence of signals also needs to be taken 
into account). The benefit of ratiometric measurement is described in the following. 
If a fluorescent molecule traverses the same path through the detection vohune many 
times (ignoring triplet-state-induced fluctuations), the number of photons detected 
from the molecule during those traversals would follow a Poisson distribution, 
characterized by a mean number of photons (this noise, which is inherent in photon 
counting experiments, is referred to as "shot noise"). This mean number of photons 
depends on the path the molecule takes through the detection volume. Taking into 
account all possible paths through the detection volume, the distribution in photon 
counts is considerably widened in comparison to shot noise. In contrast, the ratio 
between the intensity of two channels for an isolated buret is less affected since the 
mean value of this ratio does not depend on the path taken through the detection 
volume (the width of the disbibution in ratios, however, does depend on the path of 
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the molecule). Measurements using the ratio between two channels (or that consider 
joint distributions for the two channels) reduce the noise due to differing paths 
through the detection vohnne, and are thus more sensitive [1 1]. 
[0017] Dual-color cross-correlation FCS: In dual-color cross-correlation FCS, 
interactions between molecules labeled with two different colons are monitored 
using the cross-coirelation ampKtude [8. 33-35]. Significant coirelation amplitudes 
result only when a diflusing species contributes to both channels. By choosing 
different-color fluorophores that can be separated into different channels with 
minimal leakage and characterizing the backgraund, it is possible to read the 
occupancy of bomid molecules directly as the amplitude of the cross-coirelation. As 
with single-chamiel FCS, the diffusion time of the complex can be extracted. The 
occupancy and dififiosion times of the fi^e components can also be exti^cted by 
analyzing the autocorrelation of each chamiel. However, it is necessary to measure 
using a different method or to assume values for the relative brightness of the 
different species in order to extract the occupancies and diffusion times. 
[0018] Ratiometric single-molecule methods Imulti-parameter fluorescence 
detection (MFD), two-dimensional fluorescence intensity distribution analysis 
(2D-FIDA), single-pair Fluorescence Resonance Energy Transfer (spFRET)]: 
Photon burst analysis based on ratiometric methods has been developed for 
monitoring FRET, polarization anisotropy, and spectral fluctations [1 1, 36]. If the 
distance between two labeled molecules is in the 2-8 nm range. FRET cin be used to 
monitor the interaction. For example, single-pair FRET has been used to monitor 
the cleavage of DNA by a restriction enzyme in solution [12]. The same ratiometric 
data analysis can also be applied to macromolecular interactions where the 
separation between fluorophores is greater than the 20-80 A nm range for FRET, 
although this has not been done. In this case, it is necessary to excite both 
fluorophores individuaUy and perform coincidence detection. 
[0OI9J Multi-parameter fluorescence detection (MFD) and 2D-FIDA perform tasks 
similar to the ratiometric single-molecule methods, with additional capabilities. 
MFD has the additional ability to measure fluorescence lifetime [9, 37. 38], and can 
also be used to obtain the brightness infonnation available with PCH and FIDA [10]. 
Originally, the single-molecule measurements with fluorescence lifetime were 
perfonned with a single detector, although now they have been extended to multiple 
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channels. 2D.FIDA is the extension of the single-chamel FIDA method described 
above to two channels. In extracting the occupancy and brightness in each channel, 
it takes advantage of both the ratiometric and brightness infonnation. It can be used 
for samples in the low and intemiediate occupancy reghnes [39]. 
[0020] The type of infonnation that is available fiom each of the previously 
described multiple-channels method is as follows: Cross-correlation FCS extracts 
coincidence, concentration and diffusion time (and other temporal dynamics), but 
not brightness; MFD can detect coincidence and extract brightness and ratiometric 
information for multiple channels, and diffusion time, however, it can only work 
with low occupancy samples; and 2D-FIDA can detect coincidence, and extract 
brightness and ratiometric information for multiple channels. It can work with low 
and intermediate occupancy samples, but does not extract difiusion time (or other 
t^poxal dynamics). 

I0021J What is lacking in the existing mediods is a way to combine the dynamic 
information available using cross-correlation FCS with the ratiometric and 
brightness information available with MFD, 2D-FIDA, and ratiometric single- 
molecule methods, while aUowing analysis to be performed at concentrations 
corresponding to low and intermediate occupancies. 

SUMMARY OF THE TNVF.NtTnTsT 
[0022] In accordance with the present invention, a new multi-dimensional data 
analysis metiiod termed PAID is provided that aUows the extraction of information 
about temporal fluctuations (diffusion), brightness, and coincidence between 
multiple channels in a unified manner, so that aU of the characteristics of each 
species can be used for identification and for separation fiom otiier species. PAID 
also presents a convenient visual representation of the data that is usefiil because it 
focuses on photon-rich regions. 

[0023] Photon Arrival time Interval Distribution (PAID) is a new mefliod for 
monitoring both molecular and macromolecular interactions capable of independent 
and simultaneous detennination of coincidence, brightness, difEusion time, and • 
concentration of fluorophore-labeled molecules undergoing diffusion within a 
confocal detection volume. The method is based on recording the time of arrival of 
all detected photons, and then forming a multi-dhneasional histogram of photon 
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pairs, where one axis is the time interval between two photons and the second axis is 
the number of other photons detected between fliat photon pair. PATO is related to 
cross-correlation Fluorescence Correlation Spectroscopy (FCS) by a projection of 
diis histogram onto the time axis. PAID extends cross-coirelation FCS by 
measuring brightness in addition to diffusion time. Data-fitting analysis based on 
Monte Carlo simulations of diffusion paths through the detection vohime is provided 
that can be used to simultaneously extract brightness, diffusion time, and 
concentration from experimental and simulated data. The performance of PAID was 
compared to other FCS-related methods, such as FIDA, PCH, FIMDA and cross- 
correlation FCS. The statistical accuracy of the parametere extracted using PAID 
compares fevorably with the FCS-related methods, while providing additional 
information. 

(0024] The present invention provides methods for analyang/monitoring the 
properties of molecules, including marcromolecule, tiiat are labeled with 
fluorophores. The method uses a detector to detect a pluraUty of photons in a 
photon burst or other photon stream. The arrival time of each of the detected 
photons is detennined. The interval of time between various photon pairs is then 
determined to provide photon pair intervals. The number of photons that have 
arrival times witiiin the photon pair intervals is also determined. The photon pair 
intervals are tiien used in combination with the conesponding counts of intervening 
photons to analyze properties and interactions of the molecules including brightoess, 
concentration, coincidence and diffusion. The PAID mathematical model and 
resulting histograms provide useful , ways of depicting the infonnation that is 
obtained &om photon streams when analysis is based on photon-pair time intervals 
in accordance with the present invention and not fixed time intervals. 
[0025] The above discussed and many other features and attendant advantages of the 
present invention wiU become better understood by reference to flie detailed 
description when taken in conjunction with tiie accompanying drawings. 

BRIEF DESCRTPTTDN pp THE DR A WTMr.S 
I0026J FIG. 1 is a diagrammatic representation of the detection of macromolecular 
intta-actions using fluorescence fluctuation signals. In FIG. lA left panel, two 
monomers and a tetiamer of a fluorescently-labeled species traverse the laser 
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excitation profile (in green). In FIG. lA right panel, corresponding fluorescence 
bursts along with their major properties: width (red arrows; related to the size of 
species); height or intensity (bhie airows; related to the number of fluorophores on 
species); and time interval between bursts from identical species (green arrows; 
related to the occupancy of species). In FIG. IB left panel, two different molecules 
(one labeled with a red fluorophore, and the other with a yellow fluorophore) 
interact to form a complex. In FIG. IB right panel, free species emit only in one 
color, complexes emit in both colors. . Coincident detection of red and yellow 
fluorescence (orange arrow) indicates a complex. 

[0027] FIG. 2 is a schematic representation of exemplary instrumentation that is 
used for two-color cross-correlation measurements. Two laser sources (514 and 633 
nm) are introduced into an inverted- microscope; reflected by a dichroic mirror 
(DM1), and focused into a sample cell assembled using two glass coverslips 
separated by a Imm-thick silicone gasket by a 1.3 NA oil-iminereion objective. 
Fluorescence fi-om the sample is collected through tiiie objective, transmitted through 
DM1, and focused onto a pinhole (PH) by the microscope tube lens (LI). After a 
second lens (L2), the fluorescence is spht by a dichroic mirror (DM2) into red and 
yeUow channels. Red and yellow fluorescence pass through bandpass filters (BPl 
and BP2, respectively), and are detected using APDl and APD2, respectively. The 
APDs emit an electronic pulse for every photon detected, and the pulses are time- 
stamped using a counter-timer board and stored in a PC (the timed photons are 
represented in the figure as bars where the color of the bar represents whether the 
red or yeUow channel detected the photon). 

[00281 WG. 3 shows how a PAID histogram in accordance with the present 
invention is generated. FIG. 3A depicts photon streams from 3 channels (start, stop, 
and monitor). Time axis is shown in discrete clock units. A fiUed square indicates 
photon detection at that specific time position. Start photons are identified with 
different colors, stop photons with different letters, and monitor photons with 
different numbers. Each photon in the start channel is paired with each stop-channel 
photon occurring later in time. FIG. 3B depicts a PAID histogimi. x^axis, time 
interval between the start and stop photons; y-axis, number of monitor photons 
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counted in tbe time interval between the start and stop photons. Exanq,le photon . 
pair the blue start photon is paired with stop photon "d» The time interval between 
these two photons is 7 clock units, and there are 3 monitor photons (numbers 1 A 
and 3) between them. The corresponding entry into the histogram in B is coded with 
the blue background and the letter "d". 



10029] FIG. 4 depicts the generation of 2D-PAID histograms with log time-interval 
axrs and log or quasi-log monitor photon count axis. In FIG. 4A. photons detected 
m the start, monitor, and stop chamiels (filled squares) in tbe time axis (shown with 
discrete clock units. HG. 3). The contribution to the PAID histogmm for tiie fir^ 
start photon, marked with the top black arrow, is shown in HG. 4B. The time 
mterval W are shown as gray and white bands between the start channel and die 
momtorchamiel. The bins are log-spaced, with 2.4,8.16.... clock units. Themonitor 
photon count bins are shown as gmy and white bands between the monitor chamiel 
and the stop channel. These bins are.also log-spaced, with 1,2.4.8.... photon comits 
The gray and white bands below the stop photon chamiel show how the time mterval 
and monitor photon count bins combine to form the two-dimensional histogram 
bms. Each bin is labeled with the nmnber of stop photons counted wiflrin the bin. 
subsequently transferred to die histogram in FIG. 4B. 

10030] FIG. 5A shows the stmcture of the arrays with quasi-logarithmic spacing 
used for fast convolutions. FIGS. 5B and 5C iUustrate how the convolution of two 
senes of such arrays is performed. In HG. 5D, the convohition of two-step 

functions is shown. 



[00311 FIG. 6 shows the effects of changes in the characteristics of difiusing 
molecules on the single-chamiel PAID histogram. The fir^t histogram (6A) was 
formed from a single-channel simulation, with a single difiusing species 
(start=stop=monitor chamiel). x axis, time interval between the start and stop 
photons; y axis, number of monitor photon counts between tiiem. The occupancy is 
c,=0.1, the difiusion time is r°=100Ms, the brightness is ^„ =50 kHz. and the 
background count rate is =OkHz. FIGS. 6B-F illustrate how changes in the 
parameters of die diffusing species, die backgrxiund comit rate, and the sample 
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composition affect the PA© histogram. Arrows indicate changes in features (see 
text for details). In HQ. 6B. the occupancy is increased 10-fold, c.'=10xc, . fano. 
6C. the brightness per molecule is increased 10-fold. g'„=lOxg^,. InFlG.6D,the 
diffiision time of the molecule is increased 10-fold, r,'■' = lOxr^ In FIG. 6E. a 
background component is added. A^. =5 kHz . In FIG. 6F, a species with 4 times tt.e 
brightness is introduced, ,,.=200 kHz. along with a background A^, =0.8 kHz 
(values chosen so that autocorrelations corresponding to FIG. 6A and FIG 6F 
overiap). m FIG. 6G. autocorrelations, corresponding to the collapse of fte PAID 
histogmms in FIGS. 6A-F onto the time interval axis, are shown for comparison 
with PCS. 

[0032] FIG. 7 shows an exemplary PAID fit &om Table 1. In FIG. 7A. the PAID 
Hstogram for the simulation is depicted (parameters as in Fig. 6A). FIG. 7B shows 
the fitted PAID histogram. FIG. 7C shows horizontal shoes of PAID histograms 
FIG. 7D shows vertical sUces of PAID histogiBms. Slices of the simulation, black; 
slices ofthe fit, red. 

(00331 FIG. 8 shows dual-chamiel PAID histograms for sinulations containing the 
species expected in an interaction study. Symbols r and y, red aud yellow 
fluorophores. respectively, if and 7, red and yeUow detector chamaels, respectively 
FIGS. 8A, D, G, and J, are PAID histograms JtFJt, JtRR, RRY, and YYR of a mixture 
containing (c, =0.05. rf> =300 ys, ,.,=5kHz, and ,.,=45kHz). iT 

(c:,=O.O5,.f = 30OMs. 9,,=50kHz, and 9„=0kHz). and A^3^ 
(^=0.05,.3-=400ms. j3.=55kHz, and ,3r=45kHz). The background in each 
chamrel was =2 kHz . In FIGS. 8B. E. H. and K. the mixture does not 

contain the complex (^-^^ pi^sent in HGS. 8A. D. G, and J. The effect of die 
absence of ^^-^^ is seen by comparing FIG. 8A with HG. 8B; FIG. 8D with FIG 8E- 
FIG. 8G with FIG. 8H; and FIG. 8J with FIG. 8K. Differences are pointed to by 
arrows and discussed in the Detailed Description ofthe Invention. Vertical slices 
fiom the PAID histograms at time intervalr =1 ms are shown in FIG. 80. F. I, and 
L. showing tiie difference between the histograms with and without^^^'. 
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10034] FIG. 9 shows two-channel, two monitor channel PAID histograms for the 
same simulations as in Fig. 8. The red channel is assign^ as the start channel and 
the first monitor channel along tiie horizontal axis. The yellow channel is assigned 
as the stop channel and the second monitor channel along the vertical axis. Because 
the start and stop channels are assigned to different channels, these histograms 
emphasize the species that emits in both channels, the complex. HG. 9A shows that 
a two-monitor channel PAID histogram is three-dimensional, so three two- 
dimensional slices at r=100ps, r = lms, and r^lOmsare shown for each 
histogram. In FIG. 9A, both free species and complex are present In the 
r = 1 ms slice, a cartoon of each diffusing species is placed next to the contribution 
from that species. In FIG, 9B, the complex is absent. 

[0035] FIG. 10 shows dual-channel PAID histograms for samples containing 
species of labeled DNA fragments corresponding to those expected in an interaction 
study. The histograms chosen are the same as in FIGS. 8A, D, G, and J, are PAID 
histograms RYR, RRR^ JUIY, and YTR of a DNA^^''"^, DNA^^'^^®, and 
jy^^cys.wcyWB ^^^^ ^ histogmms of a 

DNA^y^'^"^ and DNA^^^'^^^ mixture. Differences in the histograms (due to the 
absence of DNA^^'^-^^^^'^^^) are pointed to by arrows. Due to the inactive 
component of Cy5, the occupancy of the dual-labeled species is significantly lower 
than the other species, lowering the peaks corresponding to the duaHabeled sspecies 
(cf.nG.8). 

[0036] PIG. 11 shows PAID histograms with SIM^RYY and fit (Table 9). FIG. 
1 1 A is the PAID histogram for the data. FIG. 1 IB is the fit to PAID histogram (1 of 
8 histograms fitted simultaneously). Horizontal slices of both are shown in FIG. 
lie, and vertical slices are shown in FIG. IID. The slices of the simulation are 
shown in black with error bars, and the slices of the fit are shown in red. 

[0037] FIG. 12 shows fluorescence-based gel images showing the major species 
present in RNAP+DNA samples. Image A shows species with TMR-based signals 
(532 excitation, 585BP60 emission); species identified as aggregates (found in the 
wells of the gel), RNAP-DNA complex and o-associated species are visible. Image 
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B shows species with Cy5-based signals (633 excitation, 650LP emission); species 
identified as RNAP-DNA complex and DNA are visible. Image C shows an overlay 
of images A (green color) and B (red color); coincidence of red and green signal is 
shown as orange. The only species showing coincidence is the RNAP-DNA 
complex. 

10038] FIG. 13 is an illustration showing flow of PAID model calculation. In the 
first step, the kernels are stretched using the scaling law in Eqs. (36) and (46). To 
account for many species the expressions k^.^^^ (A|r) are calculated using Eqs. 
(73)-(76), where the convolutions are perfoimed using the method in section 2.4.4, 
Equation (77) is used to sum over all species, and then the final expression is 
obtained by converting fi-om cumulative intensity space to photon count space 
(described at end of section 2.4.3). During a model fit, this final expression is 
compared with the data (for example by calculating x^). then the model parameter 
are adjusted. 

DETAILED DESrR JPTION OF THE JNVENTTON 
[0039] PAID in accordance with the present inverition is based on the discrete 
format of single photon-timing data (Fig. 2), which treats signals as t)hotons streams, 
rather than analog signals. Using a dual-channel microscope with avalanche 
photodiode detectors (APDs), streams of electronic pulses corresponding to single 
photons are generated, and the arrival time of each pulse is recorded, resulting in 
streams of photon-arrival times, where is the arrival time of the photon fsom 
APD^channel A, The photon stream is represented as a sum of Dirac delta functions 
(which allow representation of discrete events over die continuous variable (/): 

^.(o=i;^('-'.,> (1) 

where is the total number of photons detected in channel A \ the number of 
P * g n times ^mjn and 

^max is n = J da^{t). All sources of 

uncertainty in specifying arrival times are shorter than the >1 ^is timescale studied 
here (uncertainty arises fix)m the Heisenberg uncertainty principle (^c-2 fs), the 
APD response time (300-500 ps), and the resolution of the digital clock 
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( A/ = 12^ ns for NI 6602 counting board that times electronic pulses)). The airival 
time of each photon is recorded as an integer t„t with =[/^/a/] , (brackets 
indicate the greatest integer function, eg [3.14] =3). The arrival time is recovered 
with accuracy A/ by multiplication with the recorded integer t^,«t,,6Jl. For 

duration J . t^t will be in the range 0.1 T . where T =[r /A/J, converting Eq. 

(l)to: 

^W = |;6(t.*^> (2) 
6 is the Kionecker delta function: 6(t.t^)=l if t=t^,and 8(t,t^)=o if t*t„. In 
the following, it is assumed that the time scales of interest are much greater than the 
time resolution A<. so most of die expressions are written with the discrete time 

variable t . 

10040] Li FCS and FCCS, the relationship between a" photon stream at one 
time/,(t)and a second photon stream at a later time 7,(t+T) as a function of the 
time interval t is measured using the cross-coirelation (auto-coirelation is a 
specialization of this function), 

C.r(-)=(/.(t)/,(t + T))/(/,(t))(/,(t)) (3j 
Assuming a stationary stochastic process, the ensemble averages denoted by the 
angle brackets can be evaluated as time aveiagesrf )^ lim —£•(•.•)• By 
switching to temporal averaging over a finite experiment of duration T , Eq. (3) is 
converted into the following histogram used to estimate the cress-correlation (the hat^ 
denotes an estimator), 

^^w=T|:|;«(T.t^-t^)y';^,Ar, 

The cross-correlation histogram is formed by comparing the arrival times of each 
pair of photons t^ andt,. and adding 1 to a histogram at the bin corresponding to 
the mtervalx = t5-t^. The cross-correlation histogram can be viewed as a 
distribution of time intervals between photon arrivals (since swapping the roles of 
Is and /, is equivalent to changing the sign of x , we can specify t s 0 without loss 
of generaUty). I, and /, are defined as the "start" and "stop" photon streams 
respectively; all quantities associated with the start and stop photon streams are 
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labeled with the subscripts 5 and T, respectively. WhUe /,(t) and /^(t) are kept 
fonnally distinct, the same photon stream can be assigned to both so that 

(*) = (*) • In ftis case, Eq. (4) becomes an autocorrelation. 
[0041] To inchide the brightness infoimation lacking in the cross-coirelation, PAID 
adds a new photon-counting dimension n to the cross-correlation histogram of FCS, 
showing the contribution of each species to the cross-correlation amplitude. The 
number of photons n„ that anive between times and are counted (Eq. (4)) in 
an additional photon stream /^(t)(the "monitor" photon stream), with photon 
arrival times t„ 



t»-i 



(5) 



It should be noted that (t) is formally distinct fiom/, and I, , but can be set to be 
identical to one or both. By counting photons in each time interval as well as 
measuring its duration, the cross-correlation estimate is spread over the photon- 
counting dimension n . providing an additional history of the interval between the 
detection of the start and stop photons. In terms of ensemble averages, Ihe PAID 
function is 

(T.n) = (/, (t)/, (t)6[n,n^ (t.t + t)])/(/, (t))(/, (.)) (g) 
This function can be graphed in two dimensions (n andr), where the additional n 
axis sUces the cross-correlation into ships with different numbers of counted 
monitor photons along the vertical axis. Rather than a photon count probability 
distoibution (flie photon counting histogram of PCIWTDA). the PAID function acts 
as a photon-count distiibution of correlations. For an experiment of duration T , 
Eq. (6) can be estimated by a digital double summation histogram, 

^•') = T ||«(-.t^ -t.)5[».n, {ts„t,)]/NsN, (7) 
The PAID histogram (-»,n) describes the distribution of count rate in the stop 
channel over the number n of intervening photons in the monitor channel within a 
time intervalTsincereceiving a start photon. For each pair of photons t^ and t^.an 
event is placed in a 2D histogram, where one axis is the time interval T = t 
and the other axis is the number of monitor photons counted n =n^. Fig. 3 shows 
how a PAID histogram is formed from the start, stop, and monitor photon streams 
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usrng Eq. (7). To capture a wide dynamic range, log or quasi-log axes are used for 
the histograms. By summing over n and comparing to Eq. (4), the following is 
obtained: 

CM-tCsn^M ' (8) 
Showing that the cross-coirelation histogram is the coUapse of C^^ (t,ii) onto the 
time interval t axis. 

I0042J For a single-channel experiment, a single channel constitutes the start, stop, 
and monitor channels, [S=7^M]; the number of photons counted between each pair 
of photons is n = /-i-l for the start photon indexed by », and the stop photon 
indexed by/. This special case of the PAID histogram is written for a finite 
experiment of duration T , 

C3^(-t.n) = T i;i;5(T,t^ -t^)8(n,/-z-l) /vj (9) 
This PAID histogram is the series of waiting-time distributions to the n* photon 
(useful concept for autocorrelations that cannot be extended to cross-correlations); 
similar constructions with a single value of n were used as criteria in sifting for 
bursts [10]. By applying Eq. (8) to the single-channel result in Eq. (9), the result is 
recovered that summing the sraies of waiting-time distributions to the n* photon 
gives the autocorrelation function [40]. 

[0043] Equation (7) is extended to inchide higher order temporal correlations, or 
more monitor photon streams, by adding factors of the fonn 6(Tj.tr^ -t^) for 
temporal conelations. and fectors of the fonn 6(»,.n^^)for additional monitor 
photon streams (S„ r^. and M^: additional fonnal photon streams; and t,: 
additional photon counting and time interval variables). Higher-order temporal 
conelations are use&l for monitoring non-Maikovian dynamics [41], and multiple 
monitor photon count axes in a single histogram can take advantage of ratiometric 
measurements for better separation of species. 

[0044] Bin specification and normalization for PAID histogram: When fonning 
a PAID histogram for the large number of photons obtained in fluorescence 
fluctuation experiments, it is more meaningfiil to place events in bins than to make a 
scatter plot. In choosing the size and spacing of bhis, one needs to consider that ' 
fluorescence fluctuations occur over a large range of time scales. To cover a large 
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lange of time scales with a minimmn number of histogram bins, log or quasi-log 
time bins are commonly used in FCS. For the PAID histogram, one chooses flie 
bins for the time intervalr axis to be log-spaced, with 10 bins per decade. To use a 
log scale for the monitor photon axis is more problematic (especially at low photon 
counts) since the number of monitor photons that arrive is strictly an integer. The 
clock time resolution A/ can be chosen to be small enough to make the integer 
nature' of the discrete time interval variable x = [r/A/J negligible in the jis regime, but 
this camiot be done with the number of monitor photons counted. Unless one is 
willing to use a spacing of bins that is extremely sparse (fectors of 2.3,4...), the 
discrete spacing of the number of monitor photons will cause log bins to be 
inconsistently occupied at low n ; some bins may not even have an integer in them. 
So. a quasi-log scale is used that is adapted from the multiple tau conelation 
technique [42. 43]. The first 16 bins are evenly spaced with increments of 1, 
(n,....,n„)=(0.1,...,15); then, with each set of 8 bins, the increment is doubled 
(different increments may be used). The next 8 bins are 
°«) = (16-I7.18-19,....30-3I) with an increment of 2, foUowed by 
K.-.n„) = (32 - 35.36 - 39,...,60 - 63) with an increment of 4, etc. On the large 
scale these bins are log-spaced, while on the small scale they are linearly spaced. In 
this way, a large dynamic range of integers can be covered in a consistent manner 
with a smaU number of bins. After placing events in the histogram bins, 
normalization is necessary to obtain C^(t.»). Fh^t, the histogram is multiplied by 
T /NgNr in Eq. (7). Second, for a bin that has time interval axis limits and 
and monitor photon axis limits that include the integers through n^,^ , one 
divides by the size of the bin (x.^ -t,„,)(„,^ -„,^ +i). The value for the bin is an 
average of (t.b) over the bin limits, rather than an integral over the bin limits. 
Thisnomializes die histogram, providing C^,, (x.n). anon-ideal representation fora 
log scale: when plotting a slice of the histogram in the log scale of the monitor 
photon count axis, for a constant time interval, it is desirable that the actual area 
under the curve coirespond to the value of the correlation C^(x) . To do this, one 
approximates die photon monitor variable n as a continuous variable „, dien 
converts to a log scale using the expression f = log,.«. The relation in Eq. (8) is 
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kept valid in the new variable. Tb. sums over n are approximated as integrals over ' 
a contmuous variable „ , and converted to the log variable C : 

C5r(T) = |:c^(T.n) 

0 

= Jrf^(nlnlO)C^(x.[«j) • 

By looking at the differentials, on^finds that by changing variables ftom « to ^ the 
amplitude of the PAID histogram is changed by a fector (hio> . So, flie histogram 
bins are in addition weighted by a fector of {n+0.5)tolO. We add 0.5 because we 
consider each bin in « as covering a range between n and „+, . and the avei^ge 
over thrs range is n + 0.5. Tins only makes a difference at low n, and causes the 
« = Obua to be weighted by the fector 0.5 rather than 0. It should be emphasized 
that tbs rescaUng is made to improve the visual display of the histogmm, and in no 
way affects the statistical characteristics of the histogram. 

[0045] Efficient algorithms for constructing PAID histograms: A direct 
algorithm to form PAID histograms from photon streams compares each photon 
fiom the start chamrel with each photon in the stop chamiel, then calculates the 
number of monitor photons that fall between each start and stop photon. If there are 
N, photons in Ihe start chamiel. and photons in the stop chamiel. there ^ 
. N^Nr entries into the histogram, making the algorithm o{n^), and limiting the 
dynamic range of the PAID histogram. A more efficient algorithm exploits &e log 
spacmg of the time interval binsx. and the quasi-log spacing of the monitor photon 
comrt bins n, (Fig. 4A; filled square placed at the integer-valued time of arrival of 
each detected photon, tf; time inten-al and monitor photon count bins 
conesponding to the first start photon are denoted by the am,w at the top). The time 
mterval bins, log-spaced with integer time intervals 2.4,8.16.. . . are shown below the 
^photon stream. The monitor photon count bins, log-spaced with monitor photon 
counts 1,2.4,... are shown below the M photon stream. The time interval and 
momtor photon count bins, combined to form the 2D-histogram bins, are shown 
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below S. The number of stop photon inside each bin is shown. These values are 
transfeired to the 2D-histogram in Fig. 4B. The algorithm: (1) consider each start 
photon arrival time tf. and search for photons in T and AT that are closest to this 
time; (2) set the current time interval (r)bin to a„ = l, and the current monitor 
photon comit (n) bin to b_=l; (3) calculate the time interval at which M 
switches to the next monitor count bin b_+l; (4)if is less than the time interval 
of the (a^ +1)"' T bin, perfoim a binary search on Tto find the photon arriving just 
after r„ , to detennine how many stop photons arrive in the cuirent bin and add them 
to the (a^,b_) bin of the histogram; advance =b^ +1. and go to step 3; (5) 
otherwise, perfomi a binary search on Tto find how many stop photons arrive up to 
. the time of the (a_ +lf . bin, and add them to the (a_.b_) bin of the histogram 
and advance a_=a_+l; go to step 3, till there are no more stop photons 
o'-(«cu„.ba„)is outside of the histogram; (6) go to step 1 till there are no more start 
photons. 

[0046] This algorithm uses the feet that the start, stop, and monitor channels are 
ordered lists (each successive photon is at a later time) by performing binary 
searches. A modified search algorithm that uses increments of increasing size from 
the initial search index to bracket the desired value before performing a standard 
binary search was found to be most effective (see description of hunt in [44].) Also, 
because of the log spacing on both axes, a small number of bmary searches ca^ 
cover a large dynamic range. The algorithm is extendable to multiple monitor 
chamiels. On a 1.2 GHz Pentium 3 based PC. the algorithm shown here is able to 
form the PAID histogram of a 10 sec data set with 3.750,000 photons in 2 minutes 
and the PAID histogram of a 30 sec data set with 150.000 photons in 4.4 seconds, 
scaling nearly Unearly with the number of photons in the data set. 
[0047J Application of PAID to fluorescent species diffusing in solution: 
Fluorescent molecules diffusing in solution are excited and detected using a single- 
molecule epi-fluorescence confocal microscope (Fig. 2). The fluorescence signal is 
split into spectral regions matching the emission spectra of the fluorophores. 
Fori'-fieely diffusing, fluorescent species in solution, the species is specified by an 
index with values a=l...F; a = 0 is associated with a constant backgrcund 
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exhibiting Poisson statistics. The tightly-focused laser excitation and detection 
pinhole define the effective detection vohune . Exc(jE) is the excitation profile 
of the laser beam, and CEP{x) is the coUection efficiency function of the coUection 
opAcs [45]. V„ is defined by the detectivity, which is the product of the excitation 
and detection efficiency profilesj*(jE)sExc(iic)CEF(i): 

(x) . For simulations, a Gaussian detectivity is assumed, 
^(5)-exp[-2(x'+3;2)/fi,2_2^2/^2j^ where fi^is the l/e" width of the confocal 
volume in the and >; directions, and lis length of the volume in the zdirection 

(^eff ). For examples, a numerically calculated confocal detection volume is 

used. 

[0048] The primaiy parameters of a difiusing speciesa extracted from fluorescence 
fluctuation measurements are: (1) c,, the occupancy or average number of 
molecules mV^ ;ihe concentration [a] is obtained by dividing [a]=c^/v^ ; (2) , 
the difiusion time of a molecule across , or thfc time at which the autocorrelation 
amplitude decays by a factor of 2 (excluding other fluctuations); for a Gaussian 
detection volume with / ^ , , where is the diffusion constant; and 

0)9aA> brightness or count rate per molecule averaged over in channels. 
The average count rate for species a in channeUis \^ =c„q^^ . For background, 
the only parameter is {k^^} , the count rate in each channel. PAID histograms are 
fitted to extract all the listed parameters for the diffusing species. Other properties 
extractable from fluorescence fluctuation measurements based on extensions of 
PAID modeling, include rates for intersystem crossing, triplet-state lifetimes, 
photobleaching rates, singlet-state lifetimes, and photon antibunching [46]. 
[0049] Application of PAID to molecular interactions: The simplest molecular 
interactions are those between two species A and B, each with only one interaction 
site for the odier. For the equilibrium, A + B, AB, the concentrations of A, B, and 
AB are related by the expression K<r=[A][B]/[AB], where Kd is the dissociation 
constant [47]. Macromolecular interactions are studied by fitting PAID histograms 
and extracting occupancy, diffusion time, and brightness for A, B, and AB. Since 
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occupandcs a« dix^y proportional to concentration of species, they can be used to 
assess the affinity of the inteiBction, e.g.. the Ka can be detennined by applying 
PAID to a senes of measurements using increasing concentrations of B with respect 
to A. For mteracting species, theF diffiising species (see previous section) are A B 
and AB. . If A and B are labeled with the same fluon>phore. they are distinguish^le 
from AB because of differences in brightness (AB cames two fluorephores whereas 
A and B cany only one), and diflusion time. If A and B are labeled with 
fluorophores with distinct emissions detected by two different detection channels 
they are also distinguishable from AB because of the difference in the relative' 
number of photons detected in each chamiel (AB cames both fluorophores and is 
detected in both chamiels. whereas A and B car^ only one and are detected in only 
one channel; see Introduction). 

[00501 The dissociation constant can also be expressed in temis of fte association 
mte constants K and dissociation rate k,, K.= V Ai. An upper bound for the 
association rate can be found considering diffusion-limited encounter rates between 
naacromolecules A and B [47]. If A and B are modeled as spheres of radii and rs, 
the diffusion-limited rate constant is given by ko=4nN,iD,^Ds)ir^rs). By 
assuming that A and B have equal sizes (..=.,). and species B is present at 1 nM 
concentration, the encounter rate of any molecule of species A with molecules of the 
=?«cies B ,s k^<lQ >100-fold slower than typical difiusion times 
Conversely, if species A is present at 1 nM concentration, tiae encounter me of any 
molecule of species B with molecules of the species A is also k^^,„<io s ' To 
obtam a significant population of AB. tiae dissociation rate k, must be similar to or 
slower than the encounter mte. So. at nM concentiations. any association or 
dissociation events are exti-emely rare within the diffusion timescale. therefore, one 
only observes A. B. and AB (but not any association/dissociation events) and can 
treat the diffusing species as "static" species. 

10051] Model for PAID: The model for die PAID function will now be described 
m detail. As an introduction, a general outline and the guiding principles of the 
model are given, which are used to fit the data using non-linear least squares 
method. 



10052] The assumption that the molecules in solution difiuse independently, with 
the time between association and dissociation events much longer than ti.e difiusion 
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time. aUows the PAID function for aU molecules to be expressed as a sum of 
convolutions of the PAID function for a single molecule with the photon count 
probability distributions for many molecules. 

[0053] The expressions for single molecules are approximated by Monte Carlo 
simulation of possible diffusion paths through flie confocal detection vohune [48]. 
A scaling law is used to model changes in brightness or difiusion time (Eq. (36) and 
Eq. (46)). Depending on the accuracy desired, the confocal detection volume is set 
as an analytical Gaussian detection volmne, as a numerical approximation, or as an 
experimentally measured detection volume, allowing the direct application of the 
expected detection volume to die model. The PAID model is able to account for die 
possible diflEusion paths through that volume, in contrast to either FIDA or FIMDA 
[24. 32]. For HDA or FIMDA, only the volume density for a given brightness value 
is required. This indeed simpUfies these two models, but comes at the expense of 
not being able to model die possible diflftision paths. 

f0054] To obtain the final expressions for the PAID function through the 
combination of the single-molecule expressions, it is necessary to compute'many 
convolutions, which produces the primary bottleneck in the model calculation. 
Because of the wide temporal and dynamic ranges over which fluorescence 
fluctuations occur, logarithmic axes are desirable. For the most efficient calculation, 
the convolution method used must work in a logarithmic domain. However, pure 
Fast Fomier Transform (FFT) methods demand hnearly spaced data, which quickly 
produces huge arrays. A method is presented below in accordance wife the present 
invention that combines the use of the FFT with a quasi-logarithmic scale, making 
the model calculation practical. 

lOOSSJ Supporting theory of model for PAID: A model of the PAID function in 
accordance with the present invention is provided for several species of diffirsing 
molecules with a Poisson background in a tightly focused laser excitation volmne. 
The PAID function is expressible in terms of the photon counting probability 
distributions and the PAID function for single molecules. The photon count 
probabihty distiibution and the PAE) function for single molecules are expressed in 
terms of four path integrals, which estimate using Monte Cario simulations of 
diffirsion paths. These pafe integrals need only be calculated once; changes in the 
diffusion time and brightness parameters can be accounted for by a scaling law. 
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This means that, although the model is not expressed in closed form, it can still be 
used in a jBtting routine. 

[0056] Precise definitions of the effective detection vohnne and the brightness 
permoleculeindetectorchanneu. are needed. As a function of a molecule's 
spatial position, x=[^,y,z], neglecting intersystem ciossing to triplet states and 
assuming that the fluorescence lifetime is ze«,. the «te of photons ^^(■c)coming " 
fiom a fluorescent molecule of species a is 

\a {x) = o-„?>„ V„^Exc(^)CEF(jc). (I , J 

/.is the excitation intensity at the center o^ the confocal volume. o-„is the 
absoipdon cross-section of the fluorophore, ^„is fte quantum efficiency^of the 
fluorephore. and d„,is the detection efficiency at the center of the confocal volume 
for detectors. These can be grouped into one parameter,", , ^ 

brightness at the center of the confocal volume for detector^. Exc(3E) is the 
excitation profile of the. laser beam, and CEF(ic) is the collection efficiency 
function of the collection optics [45]. If there is more than one laser source, then 
Exc(x) may be different for each. CEF(x)in general varies as a function of 
detection wavelength. For simplicity, one uses the same Exc(-c) and CEF(3c) for 
each excitation and each detection wavelength. The detectivity is defined as the 
product of the excitation and detection efficiency profiles ^(ic).Exc(ic)CEF(ic) 
For simulations. a Gaussian detectivity is assmned 
<*(^) = expL-2(.^+y)/.^_2.V/^J, where a,is the l/e^ width of the confocal 
volume in the x and y directions, and /is length of the volume in the .direction. 
For the examples, we use the confocal detection volume for the oil immen^ion 
Objective described below. With the new definitions. Eq. (li) can be written 

[0057] THe brightness per molecule is defined as the average photon count me 
indetectorchamiel^ over the detection profile«>(x), . 
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where dV =dxdydz is an mfinifesimal volume element. For the Gaussian detection 
profile, this equation reduces to = 3^ . 

[0058] The effective detection volume is defined as in [20, 22, 49]: 

For the case of the Gaussian detection volume above, this equation 
becomes f'^ ^jr^at^l . with this definition; Eq. (12) is now written 

[0059] The occupancy c„ is the average number of molecules of species a witiiin 
the detection volume. If tiiere are N„ molecules of the diffusing species a > 0 in 
solution, and the volume of the sohition is^,, . then tiie following relationship 
with the occupancy holds, = -fs- 

10060] The density of molecules as a fimction of spatial position for species a 
« A (^) = ^ • To calculate the average count rate fiom diffusing speciesa one 

integrates the density of molecules multiplied by the intensity as a function of spatial 
position, 

^'^^=jPam^A^)dy = c^q„,. (16) 
The average count rate k,^ fiom a diffusing species a in detection channel A is tiie 
product of the occupancy c„and tiie average brightnessg„, . The total average comit 

rate in detector ^ is the sum of tiie average rate of photons of aU species and 
background 

N 

(17) 



a-1 
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Aveiage intensities are denoted by* , and instantaneous intensities are denoted hyA . • 
[0061] Spatial distribution of molecules upon detection of start photon: The 
first step in modeling &e PAID function is to obtain expressions for the spatial 
distnbution of the fluorescent molecules at the time a start photon is received. At 
the tune a start photon arrives fiom a difB^ing source a>0 (defined as time 
mtervalr = 0), the molecule that emitted the photon is inside the effective detection 
volume, meaning that the spatial probability distribution for that molecule matches 
the excitation-detection volume. This can be seen as foUows. Consider the intensity 
m the start detection chamiel of a molecule of species a as a function of spatial 
position, given by Eq.(15) while specifying the detector chamieU = 5. 

The probability that a start photon is received from the molecule at a specific spatial 
positron IS directly proportional to the intensity in the start detection chamiel of the 
molecule at that position. So, the probability distribution P„r {x) for die position of 
a molecule of species a when it emits a photon can be obtained by dividing the 
intensity 4^ {x) by its integral over all space, 

where dV=dx^dzis an infinitesimal volume element. This expression assumes 
that there is no triplet state sannation. and that ther^ is no significant diffusion 
withm the fluorescence lifetime. The superscript W means correlated, signifying 
that the subsequent photons coming fiom this molecule ai« correlated with the start 

photon. 

I0062J In contrast, the molecules of a diffusing species a>0 that did not emit the 
start photon are equally likely to be anywhere in the solution atr = 0. so the 
probability distribution for these molecules is 



V_.- 09) 



The superscript «unc» means unconelated. signifying that the photons coming fi-om 
ftese molecules are uncorrelated wift the siart photon. 

10063] Calculation of the c^oss-calculation: The cross-correlation measures the 
average rate of photons received in the stop detection chamiel T upon receiving a 
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photon in the start detection channel^, nonnalized by the average comxt rate in the 
stop detection channel. As the nmin task in calcukting the cross-coirclation, we 
obtain the average rate of receiving stop photons from a source p a time interval r 
after receiving a start photon from source a, This is the rate averaged 

over all possible initial spatial positions as well as all possible later spatial positions. 
Although this wiU allow one to calculate the cross-coirelation function, it will not 
allow the calculation of the PAID fonction without the further development of the 
theory as set forth below. To obtain the total density of molecules of a species p at 
time intervalr = 0 . upon receiving a start photon from a molecule of speciesa . one 
adds the "correlated" and "unconelated" components: 

Pas^f, (35. r = 0) = Vr {x, r = 0) 5 {a, p) + N^P /« (i, r = 0) (20) 
The second term on the right is the contribution of the molecules of species p 
that did not emit the photon. Ifa = ^ . then the molecule that emitted the photon 
gives the contribution indicated in the first term on the right, and there are N . -1 
uncoirelated molecules. SinceN^ l,onecanapproximateN^-l«N^. 
I0064I In order to obtain the cross-correlation from this mitial distiibution, the 
density of molecules at a later time can be calculated by integration againsi the 
Green's function for three-dimensional diffusion of species p [50], 

^.(-.^^) = (4.M"''exp|-(£^|. • (21) 

L ^^fi'^ J 

The resulting density of molecules as a function of the spatial variables 5 and time 
interval r is 

P<.s^p{x,z)=l dr^p^^p{x^,r = 0)g^{x,T^^), (22) 
or substituting Eqs. (18)-(20), and simpUfying, 

The total count rate on the stop detector T for photons coming from molecules of 
species p given that the startphotonwas from speciesa. *.,^^(r). is equal to the 
integral over space of the number density for molecules of species p multiplied by 
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the intensity in the stop photon stream of a molecule of species ^ as a function of 
spatial positionA^ (-f) : 

Ks-^pT (j-) = J dVp^_^^ {x, t)A^ (x). (24) 
Using Eq. (15) and Eq. (23), this expression becomes 

The first t«m is the contribution of the uncoirelated molecules of species ^ that did 
not emit the start photon. The second term is the contribution of the molecule that 
emitted the photon, ifa=^. This expression is valid for all diffusing 
speciesa./?>0. If the start photon came fix,m the background pn,cess so 

thata=0.there is no difference between the comit rate at an arbitrary time / and at 
the moment the start photon is received. So, for;ff = 0. fbe somce a of the start 
photon does not matter, 

^as-*or (^) = ^or (26) 

For^>0. Eq. (25) is still vaUd if the start photon is fiom the badcgipund 

soiu-ce,(3r = 0- 

10065] Equations (25) and (26) give us the average count mte in the stop detector 
channel T given that a photon was received in the start detector chamiel S fiom 
sourcea . These can be used to calculate the cross-con^ktion function.by summing 
*^(r)over aU somces^, as well as all possible sources a for the start photon. 
Each temi is weighted with the count rate of each source a in the start detector 
channel^, . To nonnalize the result, one divides by the product of the total count 
rate in the start chamiel k, multiplied by the total count rate in the stop chamieDli. . 
as in Eq. (3). This gives the cross-correlation fimction for channel S and channel^. 

M M 



kl. (27) 



If a Gaussian detection volume is used for the detectivity j> (if), fte standard 
fonnulas for FCS are recovered. 
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100661 While the above arguments using the density of molecules pn>vide the 
correlation function, they will not be sufficient to calculate the PAID 
fi™<^tionC^(n|r). The initial density of molecules given by Eq. (20) is still 
conect. but the information contained in the density of molecules is not sufficient to 
detennme the distribution of conelation in the monitor photon count axisn To be 
able to model the PAID function, it is necessary to know the history of the 
molecules, the possible diffusion paflis of the molecules. 

[0067] Photon count probability distribution for a single path of a single 

molecule: Before calculating the PAID function, the photon count probability 
dastrAution is obtained for the monitor detection channel given that a start photon 
was received from a molecule of species a at time intervalr = 0. A possible 
drfiusion path /of a molecule in solution is defined by its spatial position as a 
function oftimeintervalr.^,(r) = Lx,(r),;,,(r).z,(r)J. S is defined as the set of 
all possible paths / from all starting positions at time interval r = 0 to aU ending 
positions at time intervals . At time intervals = 0, the probability distribution for 
the imtial position is given by Eq. (18) for the molecule from species a that emitted 
the start photon P(x,(0)) = P;""(5,(o)) and by Eq. (19) for the molecules of 
species^ that didnot,P(j,(0)) = P7c(^^(0)^ Each of the expressions for photon 
count probability distributions in this section and the following section can be 
apphed to both correlated and uncorrelated molecules, so the superscript ' W or 
"unc" is omitted. When the combination of d.e distributions of all molecules is 
descnbed. these superscripts will be present. For species^, tire probability for the 

whole patir up to time interval r is given by the following expression (following Eq 
8 in [51]): v g q. 



= P(^/(0))exp 



(28) 



To specify tirat correlated paths are used, the superscript ' W is added. To specify 
ti^t uncorrelated paflrs ar« used, a super^ript "unc" is added instead. The 
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piobabdity for the path / is obtained by multiplying the probability density for the 
unhal position by the transition probability to every subsequent position spaced by a 
tune mtervalAr. The transition probabihty is the Green's Amotion for three- 
dimenaonal difiusion given in Eq. (21). 

[0068J The intensity in themonitor detector channel M for the molecule of species 
^ as a function of time is found using Eq.(15), 

^pM^hul^i^)]- (29) 
This expression translates, for a path/, a diffusion path into an intensity path Note 
that this expression gives the instantaneous intensity at the spatial positions along a 
drfiusron path, not the inter^sity averaged over all possible spatial positions as was 
used in the previous section. 

[0069] For a given path / of one molecule, the probability that n photons have been 
received at a time interval r in detection cham^el M is%^ The 1 indicates 

that it is for a single molecule. A differential equation-based apprx,ach is taken to 
obtam the distribution of photon counts, modified fiom [52], which contains the 
same mformation as the approach in [53J. The differ^tial equation governing the 
time evolution of 'P^^ (n |r) is: 

q; =^,m, (r) (r) 'P^ (n|r). po) 

This equation gives the rate of change in the prx>babiUty that n photons have been 
detected at a time r . The first term on the right gives the increase in this probabihty 
duetoaphotondetectedwhenn-lhadpreviouslybeenreceived. Hie mte of such 
photons is given by multiplying the total count rate A,^ (.) by the probabihty to 

have detected n-1 photons,'P,,,(n-i|.). The second tenn gives the decrease in 
the probability to have detected n photons due to a photon detected when n had 
prevrously been detected. The rate of such photons is given by mulUplying the total 
count rate Z^^(t) by the probabihty to have detected n photons, 'P^^(n|z-). The 
initial conditions are'P,,,(n = Ofr=0) = l. and 'P,^(n|r = 0) = 0 for alln^O. 
Equation (30) is an infinite series of coupled first order differential equations. By 
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Gm W - g 'P^^ (n|r) one obtains fee Poisson distribution: 
where 



(31) 



r 



C32) 



s ly paths are taken mto account the photon counting 
pn>babxlxty distribution for a molecule of species^ 

i>^(n(rj, IS a weighted 

'»V(»W-/,D(OP,(;|,).p^(„|,j ^^^^ 

i, v„.eA,„„. ,,,, ^ ^^^^ ^_ ^^^^ 

the probability for a given value of A^^ (r) , 



'^W (»(^}=/'=^A'P^ (A|r)Poi(A.n) 



(34) 
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The expression'?,^ (n|r) is the Poisson transfoim of'P^^ (A|r) . There are two 
possible spaces to work in, the photon count n -space and the cumulative intensity 
A -space. Functions in A space are denoted with a tilde. The probability density 
' ( A |r) is given by the expression 

(^K)=/,D(/)P,(/|r)^(A-A^^(r)) (35) 

This is the probabUity density, considering aU possible paths, to have a particular 
vahie of the cumulative intensity A at a time intervals . This integral is 
approximated later using Monte Carlo sampling of possible difiusion paths. (NOTE: 
this function 'P,^ (A|r) coiresponds to the probability density P (E) in [54]). The 
f^otion %^ (a\t) depends on the diffusion time and brightness of each 
species. Working in A space is advantageous because changes in flie brightness and 
diffusion time parameters for the species j3 can be taken into account by scaling in 
the appropriate dimensions. A change in g^^ corresponds to, a scaling in the A 
dimension since A^„,{t) is directly proportional to^^^. A change in 
corresponds to scaling in both the A and t dimensions. The diffusion time enters 
the model only as a product of the form D^t (see Eq. (28)). -A scaling in the time 
interval r axis accounts for a change in difiusion time in that expression. The 
scaling in A is also necessary since A^^(r) is a cumulative integral over time 
which increases proportionally with a scaling in the time interval axis. If 
%M is computed for a standard species a vidth flie diffusion time r° and 
brightaess^^^ , then for a different species fi with diffusion time t° and 
brightness ^^j^ , we have the scaling law 



So, if one calculates 'P^J^ (A(r) for the molecule that emitted the start photon, and 
'PJS?(A|r) for the molecules that did not, any differences in diffusion time and 
brightness between molecular species can be calculated using the scaling law. 
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[00711 For the background source^ = 0. assumed to be a pure Poisson source with 
intensity k,^ in the monitor detection channel, the photon count probability 
distribution is 

Po^(n|r) = Poi(A^T,n) ^7) 

In A -space, diis is 

[0072J Photon count probability distribution for all molecules in solution: The 
photon count probability distributions for single molecules is combined to obtain the 
photon count probabili(y distribution for all molecules in solution. The molecules 
are assumed to be non-interacting (at least within the diffosion time), so one can 
assume independence when the photon count distributions are combined This 
means that the photon count distribution for ±e combined source of aD molecules 
and background is the convolution of the photon count distributions for all 
molecules and background (section 2.4.1.8). The photon count distribution of 
probability for aU molecules in solution given that molecule A = l of species a 
emitted a start photon is then. 



* "P<S7 *hC^^ *)'''""'"'](n|r) (39) 



Note that the symbol for convolution * is inside the parentheses, indicating repeated 
convolutions rather than products. This is the convolution of the distributions for the 
background, the molecule that emitted the start photon, and all the other molecules 
starting with species/?=l, afl the way up to^=F. If the start photon came fiom 
the backgrounda=0. then p- is removed from the above successive 
convolutions. Fora = ^, there are N„ -1 convolutions for the uncoirelated 
molecules, since one of molecules emitted the photon and is taken into account 
by'P^S-. As shown m section 2.4.1.9. these convolutions can be performed 

equivalently m A -space and n space. To convert Eq. (39) to A -space, replace each 
n with a A . and place a tilde over each quantity? , 



m-[hu * 'P.r *n('P.r *)'^'''"-''^](A|r) (40) 
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[0073] To obtain (he photon count probability distribution for an arbitrary start 
photon, one sums over aU possible sources a of the start photon. 



P-(H-)=i;f^P^(n|r) (41) 

so] /tj ^ ' 

The weighting fector is the probabiUty that the start photon came fiorn 



source or. 



[0074] FIMDA uses a series of photon counting histograms with different time bin 
widths to extract the occupancy, difiusion time, and brightness of several diffusing 
species simultaneously. A modification of Eq. (39) can be used to model tire 
FIMDA histogram. In PAID, a photon is received at the start of each counting 
mterval. Because of this, there is a distinction between correlated and uncorrelated 
molecules. However, in FIMDA. the start of the counting interval is uncorrelated 
with the photon sequence. If the correlated photon count distribution in Eq. (39) is 
removed, one obtains a model for the FIMDA histogram. The single-molecule 
photon count distributions used in diis model are calculated as described above. 
[0075] Photon count distribution of stop channel Intensity for a single molecule 
in solution: To model the PAID function, one can combine the expressions for the 
photon count distribution wifli the cross-correlation given in the previous sections. 
The primary task is to calculate the distribution of the intensity in the stop detection 
chamiel Tover the monitor photon count variable n at a time intervalr. 
*«5^^(n|r), given that a start photon was received from source a at time 
intervalr = 0. This distribution. K,^^,{n\r), is related to the PAID function 
considering only a single species p by a constant factor: a normalization is appUed 
*° Ks^pm (n(r), similar to Eq. (27), 

[0076] For a single path / of a molecule fiom souice^ff. the monitor photon comit 
distribution of the stop chamiel intensity is calculated by multiplying the rate in the 
stop channel by the probability to have received d monitor photons, 

(n(^) = [AA^m iy) (42) 
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The superscript 1 indicates that the expression is for a single molecule. As a 
function of time intervalr. the total intensity in the stop channel is^^(r) . 
\^ {^\r) is how this intensity is on average divided up among the different values 
of the monitor photon counts n . 

I0077J When all possible paths are taken into account, the monitor photon count 
distribution of the stop channel intensity for a single molecule of species^. 
'W(nk). is a weighted average over all paths of the distribution for a single 

As in Eq. (35). P^ (/(r) is given by Eq. (28). and 5 is the set of paths fiom all initial 
positions to aU final positions after a time interval . . Eq. (43) can be rewritten as an 
mtegration of die Ppisson distribution against the contribution to the stop chamiel 
mtensity for a given value of cumulative intensity A , 

00 

'kp^ {^\^)=ld\%„ (A|r)Poi(A,n) (44) 

where 

%»a (A|r)=/^D (/)P, {l\r)A^ {z)S{A-A„^ (r)) (45) 

The difference between this equation and Eq. (35) is the additional factor of the 
intensity in the stop channel of fee path/,^^(r). Because of this fector. Eq. (45) 
gives the cumulative intensity distribution of the stop channel intensity, rather than 
the cumulative intensity distribution of probabiUty. 

10078] The scaling law in Eq. (36) also applies here, with one additional fector. the 
brightness in the stop detector chamiel^^. if %^{^A\r) is computed for a 
standard species a with the diffusion time r« and brightness in the monitor 
chamiel and brightness in the stop chamiel, dien for a different species widi 
diffusion timer|. brightness q^^ in the monitor chamiel and brightness in the 
stop channel, we have the following scaling law 
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(46) 



As w,th the cunndative intensity probability distribution in Eq. (35). one evaluation 
of the cumulative intensity distribution of the stop channel intensity is made for the 
molecule that emitted the start photon, '^^ (a|.). and one for the molecules that 
did-t.'^- (A|r). Any changes in the pa^nete,..;. ,,,.and can be taken 

into account by using the scaling law. 

10079, For the background source^ = 0 . assumed to be a pure Poisson source with 
mtensxty k,, in the monitor detection chamael and intensity ^ m the stop 
detection cham^el. the photon count distribution of stop channel intensity is 

W(n|r)=VPoi(Ai,j^T.n) ^^^^ 

In A -space this becomes, 

[0080, PAm function for all molecules in solution: To account for all of the 
other molecules in solution. 'A^^Cnlrjis convolved ^th the photon count 
probability distribution for all other molecules and background. If the molecule in 
question emitted the start photon, then 



The fu.t factor is the background photon count p„,babili.y distribution. The second 
factor IS the monitor photon count distribution of stop chamiel intensity for the 
molecule that emitted the start photon. n.e r.st of the fectors come fiom the 
molecules of all the species that did not emit the photon. These are successive 
convolutions, not products, as indicated by the star * inside the parentheses 
Fora =y , there are N„ -i convolutions for the unconelated molecules, since one of 
molecules emitted the photon and is taken into account bv'*«" 
[0081, If the molecule in question did not emit the start photon, then 
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(N/r -Hoc,fi))^P^ * -p- • »]^(.p«» j^^i^^ (50) 

The pre-fector is how many molecules of species there are in solution. The firet 
convolved fector is the background photon count probability distribution. The 
second convolved fector is the monitor photon count piobabiUty distribution for the 
molecule that emitted the start photon. The third convolved factor is the monitor 
photon count distribution of stop channel intensity for the molecule of species ^ 
that did not emit the start photon. The rest of the convolved fectois come from the 
molecules of all the species that did not emit the photon. Fora = y . there is one less 
convolution for the uncorrelaled molecules, since one of molecules emitted the 
photon and is taken into accomit by'p- . Forfi^y, there is one less convolution 
for the uncoirelated molecules, since one of molecules is the one whose intensity is 
being calculated and is taken into account with'A^ 
[0082] For the background. 



(51) 



The first factor is the monitor photon count distribution of stop channel intensity for 
the background. The second fector is the monitor photon count probability 
distribution for the molecule that emitted the start photon. The rest of the fectors . 
come from the molecules of all the species that did not emit the photon. Fora , 
there is one less convolution for the miconelated molecules, since one of molecules 
emitted the photon and is taken into account by'P^ . 

[00831 The contributions to A^^^(n|r) of all molecules is now summed, 
obtaining, 

*a..;^(n|r) = A-^(„|r)+A-_^^(n(r) (52) 

The first term on the right is for any contribution that is correlated with the start 
photon, and the second term is for the contribution that is unconelated with the start 
photon. To obtain the final expression for the PAID function, we sum over all 
possible sources a of the start photon, weighted by the intensity of each somce in 
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the start channel, just as with the cioss-conclation. To nonnalize the result, one 
divides by the product of the total count rate in the start channel multiplied by 
the total count rate in the stop channelA^ . as in Eq.(27). 



U if 



I]Z*a5*<.s-»/Wtf(n|r) 
C^inlr)-^^^ 

Remember, to convert any of the n -space expressions to A -space, replace each n 
with a A . and place a tilde over each quantity P or* . Uris is the final expression for 
the PAID function. It is calculated by adjusting the distributions 'P™ (A |r). 
'^^(Ak). 'Pr(A|.) and ^kZ.{A\r) for the parameter, of each species^, and 
usirig Eqs. (49)-(53). A more detailed description ^vill be prx>vided below tfaat 
demonstrates how to calculate the PAID function. 

[0084] Combination of sources for probability distribution: As an example, 
assume one is given two photon count probability distributions. P, (n,)andP, (n,) . 
and one wants to determine tte photon count probability distribution fL 
combinedsomceP(n),wheren=n.+n,. This can be determined by 

P(°) = SP.K)B,(n-n,|ii,) (54) 

To obtain the probability to count „ photons fiom the combined source, one smns 
over aU possible values of n.. For each valuen. , the value of n, is restricted so that 
the value of the combined total counts is conect,n, =n-n, . The probabihty to 

have n, counts from somx« 1 isP,(i,,). and the probability to have n,=n-n. 

comitsinsouiceZgiventhatthereweren. counts fiom source 1 is P, (n-n, |n,). K 

is assmned that sources 1 and 2 are independent; so thatP, (n -n. | „.)= p, , 

and Eq. (54) reduces to a convolution. 

I'(°) = S;P. (n.)P,(n-n.).(P. *P,)(„) ^^^^ 

To combine more than two independent sourx^s. one simply performs successive 
convolutions, P (n) = (P, *P,*P^ * • )(»)• 



wo 2004/011903 



PCTAJS2003/023252 



-39- 



(0085J Equivalence of convolutions in n-space and A-space: One property 
tbat needs to be established is the equivalence of performing convolutions in n 
space and Aspace. Consider three distributions f(n|r). g(n|r), andh(n|r). such 
that h is the convohition of f andg : 

l^(»k) = Zf(»10g("-»'k)-(f*g)(n|r) (56) 
One can write each of these distributions in n as a Poisson transformation. 

r(ii|r) = /f(A(T)Poi(A,n)rfA 

0 

00 

g(n|r) = /g(A|r)Poi(A,n)rfA (57, 

CO 

h(n|r)=Jh(A|r)Poi(A.n)rfA . 

0 

By substituting the Poisson transformation expressions for f (n.r) and g(n.r) into 
Eq. (56), one gets 

h(n|.) = |^A,f(A,|.)/rfA,g(Ajr)Jpoi(A,n-)Poi(A..n-n-) 

= JrfA,f(A,|r)/rfA,g(Ajr)Poi(A,+A,.n) (53) 
= ^^[pA'f(A'|r)g(A-A'(r)jpoi(A.n) 

Smceitis knownthath(n,r)=JrfAh(A.r)Poi(A.n).onefinds 

0 

h(A|r)=|rfA'f(A'|x)g(A-A'|r) = (f*g)(A|r) (59) 

This shows that convolutions can equivalently be perfonned either in n -space or A - 
space. 

IOO86J Combinations of sources for the monitor photon count distribution of 
stop channel intensity: Tie way to combine the monitor photon count distributions 
of stop chamiel intensity for multiple sources will now be demostrated. This 
distribution, k^^^ (n |r) . is related to the PAID function considering only a single 
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species p by a constant factor a nonnalization is appbed to A:^^^(n|r). similar 
toEq.(27) •c^(„|.)=:^£^£^Hl) consider two intensity paths and # 

fiom independent sources 1 and 2, respectively. The intensity of each path as a 
function of time interval is given for each detector channels. and;ij,(r). 
Each source may be a single molecule, or more than one molecule. For each path, 
the monitor photon count distribution of the stop channel intensity is calculated by 
multiplying the rate in the stop chamiel by the probability to have received n 
monitor photons, 

C,(n|r) = Pi,(n|r)4,(r) 

For the path r, of source 1, the total intensity in fee stop channel is .l^^(r) as a 

function of time intervals. *i^(n|r) gives how much of that intensity on average 

is detected foragivenvalueofmonitorphotoncountsn. >t^^(n|r) is interpreted 
similarly. 

[0087J. When aU possible paths are taken into account, the monitor photon count 
distributions of the stop chamiel intensity. C(n|r)andA^(n|.). are weighted 
averagesoverallpathsof ^i^,(n|.) andA^,(n|r). respectively. 

*-(°W=t.D('7)P'(7|r)C(n|r) 

P'(7|r)is the probability for a given intensity path^, and is the set of all 
intensity patiis fiom up to a time intervalr. P*(^|r) and 5^ ai« defined 
analogously. The monitor photon count probability distributions Pi(n|r) and 
K (n I'-) are calculated by similar weighted averages, 

P^(»k) = L.D('7)P'('7W,(nk) 

Pv(nk)=i,D(^)P^(^|.)Pi^(„|,) (62) 
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r0088J When one desires to determine the monitor photon count distributions of 
stop channel intensity for the combined source. 

As before, the monitor photon count distribution of the stop channel intensity is 
calculated by multiplying the xate in the stop chamael by the probability to have 
received n monitor photons, 

W(»k) = P«^(n|r)^^(r) ^g^j 
[0089] When all possible paths are taken into account, the monitor photon count 
distributions of the stop channel intensity for the combined somce, (n|r) , is 

Now. since sources 1 and 2 are independent, the monitor photon count probability 
dxstnbution for the combined somce is the convolution of the distributions for the 
individual sources. 

P«^(»W = (Pi.,*P^,)(n|.) (6,^ 
Using this along with Eq. (63), Eq. (64) is rewritten. 

By substihiting this into Eq. (65) and grouping tenns fiom the same source, one gets 

= D (7)P ' (n •|r)J 

Now, using Eqs. (61) and (62), one gets 

V(°k) = (A];^»P^)(n|r)+(*^*P^)(„|^) (ggj 

This equation convolves the monitor photon comit distribution of stop chamrel 
mtensrty for each soun.e with the photon count distribution of the other, and then 
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adds the results. This equation alone with En u 

^uauuu dxung wim jiq. (66) can be appbed successively to 

combine many sources. 

This expx^sion is used in this description to combine the monitor photon count 
distnbutions of stop channel intensity for aU of the molecules in solution 
[00901 Evaluation of kernels for model: Tt. distributions'P-(A|r). 
'^^m. 'Pr(A|0 and ^k^(A\r) for the standard species ^ ax. evahiated 
using Monte Carlo generation of possible diffusion paths /. a modification of the 
simulations described. Tte standard species has a diffusion time^ =1. a brightness 
= 1 in the monitor channel and brightness =1 in the stop chamiel (each of 
d.ese has arbitrao' units.) As in the previous section, a diffusion path / is simulated 
by a series of three-dimensional r^dom distance steps, with mean ^ = 0 and 
standard deviation a = 72^, ^fa^re D is the diffusion constant and A. = 10-.'' 
is the time interval step. The initial positions a« drawn from a unifonn distiibutiol 

across the simulation box P (x fO^^ - ^ r^*u . . 

v^' ^")) - y «>r the uncorrelated molecules, and from 

box 

the detection profile P (^,(0)) = P^~^(,,(o))for the conelated molecules. The 
uncorrelated distributions are formed from simulations restricted to a finite box of 
size with periodic boundary conditions. Since the probability density for the 
initial uncorrelated spatial position is a constant with respect to spatial position, the 
srze of the simulation box must be restiicted to have a significant number of 
drf&sron-driven crossings of the detection volume. TTae periodic boundary 
conations are necessa^r to avoid a dr^ in tire number of ti.e molecules from the 
smxulatioa box. For the correlated distribution, there is no restiiction on tiae 
smxulation box size. The probability density for the initial correlated spatial position 
rs restiicted to tire detection volume, and the exit of tire molecule fiom tire region of 
the detection volume causes ti,e average number of molecules per detection volume 
to return its uncorrelated vahie. For tire speciesa. the joint probability for the 
whole patir up to time interval ris given by tire foUowing expression (Eq. (28) 
without the limit): ^ m v 
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10091] The distributions are evaluated at a series of time intervals r that is 
logarithmically spaced overlO decades fiom r = 10-.f to r=10^r- with 100 bins 
per decade. The cumulative intensity A bins are logarithmically spaced over 10 
decades fiom A = 10-g,,r- to A = 10^^,,.^- with 10 bins per decade. At each 
time interval r in the logariftmically spaced series, the cumulative intensity 
KA'^)^lKM,{r')AT' and the intensity in the stop chamiel (r) ar« evaluated. 
An entn. of 1 is entered at the corresponding (r. A) bin in the cumulative intensity 
pn>bability distribution. ^f^{A\.)or^f:^^i^A\r). Also, an entry is added 

to the corresponding (.,A) bin in the cuxnulative intensity distribution of stop 
chamael intensity. '>fc.'S^(A|r)or'^„":^(A|r). After simulating a total of 10* paths in 
the correlated case, and 10^ paths in the uncorrelated case, the estimated 
distnbutions are divided by the number of paths used. In the uncorrelated case 
there are actually only 10< diffusion paths simulated. Each diffusion path is used 
withlO' evenly spaced starting points to give a total of 10' effective paths. We can 
nse many different starting points in this case because the molecules are 
uncorrelated: there are no special properties of the initial position. 
[0092] Implementation of the model and fitting routine: The following is a 
summary of how to calculate the PAID function for a given set of parameters for 
Fdiffusing species and background: the difiusion time.-, the occupancy^, . the 
brighti^ess per molecule^^, . and the background intensily^„, . For each sp Jes^ 
one calculates the single molecule cumulative intensity distributions of probability 
and stop chamiel intensity for both ti.e correlated and uncorrelated initial positions 
Eq. (36) is used to calculate 'P,-^(a|.) and •p^;?(A|r) and one uses Eq. (46) to 
calculate '^,^(A|r) and '^-,(a|.) from the distiibutions estimated for the 
standard set of parameters, m this way. the difiusion time of each speciesr-. 
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toightnessinthenioiritcM-charmel^^^, and brightness in the stop cte^ are all 
accounted for. 

[0093] m building up the PAID function (n|.) from these initial distributions, 
many convolutions are necessary. In order to perfonn them quickly, we have 
developed an efficient algorithm to compute convolutions in a quasi-logarithmic 
scale (see next section). "P,"^(a|.). 'P,^(A|.) and '^|^(a|.) are 

rebimied into the quasi-logarithmic scale described below, and the convolutions are 
performed in A - space. 

[0094] Now. Eqs. (49H52) are used to calculate the monitor photon count 
drstnbutions of stop chamael intensity for all of the categories of molecules For 
each difiusing species;?>0. there is a distribution for the correlated and 
uncorrelated molecules. 

[0095] nre simulations used to form the distributions 'P-(a|.) and 'k^{K\r) 
for correlated initial positions are not restricted to a finite simulation box. Initially 
the molecules are inside the detection volume, but the molecules are allowed to' 
drfiuse out of the detection volume without restriction as the simulation time passes 
Because the size of the simulation box is um«stricted. the difference between 

-6(a.r) is ignored: in Eqs. (49)-(52). the expression 6(a.r) is dropped. ^ 
[0096] The simulations used to form the distiibutions ^f^W and ^k^(A\r) 
for uncorrelated initial positions are restiicted to a box with a volumeF^ . which is 
larger than the effective detection volume: K^, V^. The box has harmonic 
boundao. conditions, so that a molecule that comes out of one side reenters the other 
side. The fundamental concentiation for the distributions formed fiom these 
simulations is^ . Successive convolution of the distribution ^f^ (A)r) models 

an increasing concentiation. but only in discrete steps. If the number of unconelated 
molecules of species inside this box isA^^^, tiren one models the uncorrelated 
photon count probability distiibution using('P-(A|r)*f- . To have the proper 
occupancy, then the following relation should be satisfied. 
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[0097] Because the volume of the simulation box used for the distributions for 
uncorrelated molecules is much smaUer than the volume of solution the 
distributions •p-(A|r) and •*,^(a|.) account for a much larger concentration 

thandescribedpreviously. J- rather 4an J- . To obtain fl.e proper occupancy; the 
number of molecules N, in Eqs. (49)-(52) is replaced by the number of molecules 
in the simulation box,Ar^. Using the small box size for the uncorrelated 
molecules and the.large box size for the conelated molecules. Eq. (49). which is 
for correlated molecules, is rewritten in A -space, 

''^-*/>Mr(A\T) = 8ia,fi)^f,^ * *fl('f^ *r'^](A|r) {73) 
Equation (50), which is for the micorrelated molecules, is rewritten 

Equation (51), which is for the background, is rewritten 

Finally, Eq. (52). which combines these expressions, is rewritten 

*«-.^(A|z-) = *-,^(A|r)+A-^^^(A(r) (75) 
If is not an integer, then the distributions ai^ calculated with the closest 

integer nmnber of molecules in the simulation hcN^^ , adjusting the brightness 
9^ to satisfy the relation g',,Ni,^ =q,,N^. This keeps the total count rate 
from the species /3 constant, while using the closest integer for the number of 
molecules in the simulation box. Note that the adjusted value for the brightness 

is used only for the expressions for the uncorrected molecules, not the correlated 

molecules. 
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[0098] Converting Eq. (53) to A-space. one gets for .be cumulative intensity 
distribution of correlation. 



M M 



'-'^^^''''^ — - 

This is then converted to n -space to get the PAID function, 

CO 

Qw(n|r)=JrfAC^(A|r)Poi(A,n) 

[0099J The kernels used as the basis for the model have finite sized bins in A , which 
are indexed by b and have the raBge[At,.AL). Because of the finite bin size, 

what is really calculated is the average over a ^'n ^^ "^^'^sru j^'lr) 

A* -A*. • " " 

assumed that^ the amplitude is constant across the bins, so that 

' ^ AL^Xl f^"-^ ^°fl>erange[Al,Al). 

[0100] The expression for C^(n|.) as calculated until now gives the 
instantaneous rates at a particularr . However, the time interval bins for the photon 
counting data have finite extent, and so (n|r) must be averaged over Ae range 
ofthe time bin. To calculate C^(A\r) within the time interval andr^.one 

inteipolates between C--,fA|r . ^ andC (Air \ n ( k\ \ • 

^\ N^nn. ; anac^ ) - (A|7^) is projected 

forward in time fromr,, , noting that the limits of a bin in cumulative intensity 

scales with time interval. A', MsA'^ — ptiHA' (^\ ^ 

» >^«a, \J) - ^rnrn and (r) = A^ . The value 

mtn 

interpolated forward fiom isC^ (A|r) ^i^^^^^^^ll^^ ^^^^^ ^ 



chosen so that A*^ <A-^i!!ss.<. a* <:;^;i„.i ^i. , 

- A ^ < A^ . Smularly, the value mterpolated backward 
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sm{^\T)~— , where b is chosen so that 

(A^ - j 

T" 

max 

A^ ^^^< AL . The intenjolated values are averaged over the rectangle 

bounded by r^, A^.andA^. Note that more than one bin may contribute 

to the averaging. 

[0101] The final integration over Asho,vn in Eq. (78) is implemented as a matrix 
multiplication. Tte app„,ximation thai the vahae of (a|.) is constant over a 
binwidilimitsA^,and A^, leads to the approximation. 

= (79) 
Hae integral on the right can be expressed in tenns of the incomplete gamma 
function;' = -i-L-f-rf^ 

JaL ^°H^°)''A = r(n + l.AL)-r(ii+l,Ati„) 

The bin spacing is kept fixed, so the integrals in Eq.(80) need to be perfonned only 
once to create the matrix. 

f0l02J Efficient calculation of convolutions in a quasi-logarithmic scale: 
Because of the large mrmber of convolutions necessary to calculate the PAID 
histogram, it is necessary to have an eflScient algorithm for convolutions. The 
standard method is to Fast Fourier Transform (FFT) the data, multiply the data in the 
complex transform space, and then FFT back to real space (See for example Ch 13 
of [44]). The FFT, however, requires evenly spaced bins, which can produce huge 
aixays if one wants a large dynamic range. A convolution medrod inspired by the 
multrple-tau correlation technique [42, 43] that uses the FFT for small convolutions 
IS preferred. For each vector (a andb) to be convolved, one produces a series of 
Inrearly spaced arrays of lengdr 64, where the resolution of each is reduced by a 
factor of 2 compared to its predecessor. The stmoture of these vector, is shown in 
Frg. 5A. m the figure, dre arrays are of length 8 (rather than 64) for each factor of 2 



(80) 
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in resolution. The anay with the smallest spacing is shown at the left. In Fig 5A 
the first array has a spacing of AA =0.1 . n,e second anay has double the spacing.' 
AA - 0.2 . The first four elements axe each made up of a smn of two elements fiom 
the first array. The third amy has a spacing that is doubled again. AA=04 
Again, the first four elements are each made of a sum of two elements fiom the 
second anay. Continuing this structure to successively lower resolution, one obtains 
a senes of linearly spaced anays that can efficiently span a large dynamic range 
I0103J m Figs. 5B and 5C. we illustrate how the convolution of two series of such 
arrays is performed. One anay is labeled with lower-case letters, and the other is 
labeled with upper-case letters. To calculate the convolution of two vectoxs a*b 
the arrays with the smaHest spacing (AA = 0.1) are convolved with each otter first' 
as shown in Fig. 5B. The FFT procedure described above is used. whe,. the arrays 
are zero-padded to twice the original length. The lower half of the array becomes 
the new array for the smallest spacing. TTae whole array is rebimred to dre next 
largest spacing, and serves as a contribution to the convolution with spacing 
AA = 0.2 . ^ a 



[0104] In Fig. 5C, it is shown how the convolution for the larger spacing AA = 0 2 
IS perfonned. Because part of the convolution with that spacing has ah^dy been 
performed, one needs to exclude that contribution in subsequent calculations This 
IS done by setting (or "clipping") the lower half of one sub-array (lo^er-case) to 0 
and convolving it with the other sub-array (upper-case). The clipped version fix,m 
a -convolved with the full versionfi.mb.and vice versa. By adding these results 
wrth the array obtained by r^bimring the result with spacing AA = 0.1 . one obtains 
flie final result for AA=0.2. TTris allows the use of infomration fix>m higher 
resolutron arrays to contribute to Are convohitions with lower resolution arrays 
[0105] In the quasi-logarithmically spaced arrays, the function to be convolved can 
be approximated by a series of steps; over the range of each bin. the fimction is 
assumed to be constant. When two such steps ar« convolved, one gets a triangle 
shown in Fig. 5D. The convolution of the two steps leaks into the next bin nre 
discrete convolutions shown in Figs. 5B and 5C do not account for this because drey 
assume that dre function is nonzero only at discrete values. In order to approximate 
the continuous fimctrons widr the method presented, dre convolution is shifted one- 
half bm up after each FFT-based convolution in Figs. 5B and 5C. 
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[0106] In Eqs (73)-(75). there are successive convolutions perfonned on the same 
anay which have the fonnCp- . in order to compute successive 
convohitiom on the original distribution'?-, we" can use a trick to prevent 
perfonning N^^ convohitions directiy. The number is expressed in binary. 
By recursively convolving the original amy. one obtains a series of vectors for each 
convolved power of 2: 1.2,4.8,... By convolving only those convolved powers of 2 
in the binary representation ofiVj^, one obtains (' •j'^ . This allows the 
number of convolutions perfoimed to increase logarithmically with 



taw 



Examples of PAID appUcations 

10107] In the foUowing examples of the present invention we use single- and dual- 
chamiel simulations and experiments to illustrate some of the appUcations of the use 
of PAID. Labeled DNA constructs as well as an RNA Polymerase-DNA interaction 
are used to show that PAID obtains expected results. These examples axe not 
limiting; many other applications are possible, including those with multiple 
channels, and more monitor channels. 

Materials and methods 

[0108] Preparation of DNA: We used fluoiescentiy-labeled DNA fiagments as 
model systems for exploring the capabihties of PAID. The fluorophores used were 
Cy3B 0^560 nm.A«.~580 nm). Cy3 (a«~550 nm,>^~570 nm) and Cy5 
(X«,~650 nm.A«„~670 nm). Six DNA fiagments were synthesized (the 
nomenclature used in the superscript report the fluorophore, the DNA position, and 
die DNA strand where the fluorophore was inhoduced; T: top strand, B: bottom 
strand): (1) DNA^y^^-'^, (2) DNA^^^''^''^^^-"^ (3) DNA'^^-'T, (4) DNA^^-^^ (5) 
DNA^'.>^«V3.65B DNA-y5.6» ^NA fiagments were prq,ared using stankrd 

PCR protocols [55] witii one or two 5'-labeled DNA primers, followed by 
purififcation using non-denaturing gel electrophoresis. The sequence had 65 base 
pairs. Mixtures of DNA^^-'^ and dNA*=^^»-'^'«^^».«b ^^^^ ^^^^^^ 
channel applications; and mixhires of DNA*^''"^, DNA^^-"^ and dNA^5 '™^3.«5b 
for dual-channel appUcations. In aU cases, the intramolecular separation between 
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tbe two fluorophoxes is lai^ge. (65 bp. -240 A) to preclude FRET between the 
fluorophores. Data wexe acquired for 5 min usirrg 30 pM-1 rM DNA in 20 mM 
HEPES-NaOH (pH 7). 50 mMNaCl. 5«/o gfycerol. and 1 mM men^aptoelbylannne 
(Fluka^ Milwaukee. WI). We added O.P/o BSA (Panvera. Madison. WI) to the 

and DNA^«-'T samples to reduce adsorption onto surfeces He 
concentration of DNAC^^-- was detemuned using UV-Vis spectrophotometcy and 
fluorescence spectrophotomeby to calibrate concentrations detennined fiom FCS 
measurements. Higher concentration (10-25 nM) samples were prepared for all " 
DNA fiagments. and occupancies and diffusion times were extracted using FCS as a 
basis for dilutions. 

[01091 Preparation of RNAP and KNAP-DNA Complex: E.cMia coH 
RNM> core was purchased fix,m Epic«rtre (Milwaukee. WI). and RNAP o subunit 

(c 0 was purified and labeled at amino acid Cys569 using tetramethylrhodamine 
(TMR) as described [56]. RNAP holoenzyme and RNAP-DNA complexes were 
formed essentiaUy as described [56]. using DNA^^-«B for the formation of the 
complex; the complex was diluted to I rM nominal concentration for the PAID 
measurements. The large distance (»100 A) between TMR and Cy5 in the RNAP- 
DNA open complex precludes FRET between the fluon>phores [57] After 
formation of Are RNAP-DNA complex, the sample was loaded in non-denaturing 
5% polyacrylamide gels and was electrophoresed at 10 V/cm for 1 hour the 
resulting gels were imaged using an x-y fluorescence imager (Molecular imager FX 
Brorad. Hercules. CA) equipped with 532-mn and 633-m„ excitation lasers (for 
excitation of "yellow" and "red" fluorophores respectively), and 585BP60 and 
640LP emission filters (for detecting "yellow" and "red" chamiel emissions 
respectively). ' 

[0110] Confocal fluorescence microscopy: The instrumentation used is similar to 
that described in difiusion-FRET studies [(see Fig. 2) and [11. 12]]. For single- 
chamrel experiments, the 532 mn line fiom a solid-state pumped Nd:YAG laser • 
(GCL-IOO-S. Crystalaser, Reno, NV; 100 kW/cm^ excites Cy3 and Cy3B) was 
mtroduced using fiber optics. An excitation dichroic mirror ^DMl) was used to 
reflect the laser light, while transmitting fluorescence emission (400-535-635 
TBDR. Omega Optical. BratUeboro. VT). For the dual-chamiel experiments (DNA 
fiagments and RNA Polymerase), two laser beams. 514 mn Ar* (543-A-A02 
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Melles-Griot, Carlsbad. CA; 200 kWW; excites Cy3 and TMR. and to a much 
lesser degree. Cy5). and.633 nn. HeNc (05.LHP.171. Melles.Griot.Carlsbad.CA- 66 
kW/cm^; excites Cy5) were nsed. A different dichioic minor (DMJ; 390-510-630 
TBDR, Omega optical. Biattleboro. VT) was used to reflect the two laser lines 
while transmitting the emissions of TMR. Cy3. and Cy5. For both experiments, the 
laser excitation is focused 20 ^m inside the sototion by a lOOx 1.3 NA Zeiss 
Neofluar oil-immersion objective. Fluorescence fiom the detection volume is 
focused by the microscope tube lens (LJ) on a 100 pm pinhole (PN) and split with a 
dachroic minor (DM2; 630 DMLP. Omega Optical) into two detection chamiels- the 
Cy5 chamiel [filtered using a 650LP filter (BPI), Omega Optical], and the 
TMR/Cy3/Cy3B channel [filtered using a 580DF60 filter (BP2), Omega Optical] 
Sdicon avalanche photodiodes (APDI and APD2; SPCM-AQR-14. PerkinEhner 
Vaudreuil. QB. Canada) detect fluorescence photons, photon-associated e'lectronic 
pulses are timed by a counter-timer board (PCr.6602. National Instruments. Austin. 
TX), and stored in a PC. 

[0111] Simulation of translational diffusion of molecules, and photon emission 

and detection: Simulations for translational diflusion. photon emission and 
detection of molecules were perfonned simUarly to [58]. A C3aussian detection 
volume with «, = 0.35 nm. / = 1.75 ^un is placed at the center of a SD-simul^tion box 
with size =3.5x3.5x17.5 nm'. assuming periodic boundary conditions (a 
molecule that leaves V^^ reappears at the opposite side with the same lateral 
position); a fixed number of molecules is placed inside the box. Diffiision in and out 
of the detection volume is simulated by a series of steps of A.=i^s, a timescale 
short enough to ignore the effects of diffusion ^ lOO for our systems) on the 
excitation rate of the molecules within one step. At each time step, the distance step 
for each dmiension (x.y.2) is determined by a pseudorandom number generated with 
a Gaussian distribution [44], with mean ^ = 0 and standard deviation <t=V25a7. 
where D is the difiusion constant. The diffusion-step-generation distribution is 
taken direcUy fiom the Green's fimction for SDnlifiusion [50]. For each difiusion 
step, a series of pseudo-random numbers is generated with an exponential 
distnbution with a decay rate X that depends on the excitation rate for the 
molecule's position until time A. is passed, thereby generating a series of arrival 
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times for emitted photons; the fluorescence lifetime is assumed to be zero, thus 
ignoring saturation and other photophysical effects. Finally, a second 
pseudorandom number is generated to detemiine if the photon is detected, and, if it 
is. which channel detects the photon. The time of arrival and detection channel of 
each detected photon is saved. The resulting photon sequences can be subjected to 
any of the possible data reduction and analysis methods. The idealized conditions - 
described above (Gaussian detection volume, neglected fluorescence lifetime) are 
used for the simulations since they are well-modeled by each of the data analysis 
methods. 

10112] Fitting routine: The Levenberg-Marquardt nonlinear least squares fitting 
procedure [59] is used to extract the parametera &om the data, calculating the 
necessary partial derivatives numerically. The model used is described in detail 
above. To estimate the standard deviation of enora. we use ten independent 
instances of the histogram to calculate an estimate of the statistical errors of each 
bin. The statistical errors are used as weights in the fitting routine (this method is 
used for PCS in [58]. If a bin is nonzero in fewer than ten instances, it is excluded 
&om the fit. Without this restriction, we found that only a few points dominate the 
value for ;jr' (the merit function is not ideal for sparwly populated bms in 
histograms). 

[0113] Calculation of detection volume and diffusion parameters: A 
numericaUy-approximated, non-Gaussian volume is used for the analysis to account 
for our experimental conditions. The laser excitation profile for an oil-immersion 
objective focused 20 ^m inside the aqueous solution was calculated using Monte 
Carlo integration of plane-wave contributions using expressions fiom [60, 61], 
simplified to scalar diffraction, and accounting for the water-glass diekctric 
interface. The effect of the pinhole in the detection path (100 ^un) was calculated 
using geometric optics [45], accounting for the water-glass dielectric surface. ITie 
peaks of the excitation profile and detection profile were translated with respect to 
each other to assure good overiap in the z direction (adjustment of the pinhole or the 
APDs along the emission light path to maximize the signal would produce similar 
results). By multiplying the two profiles, we obtain the detection vohmie which 
forms the basis for calculating the kernels used for fitting the experimental data. 
The size of the calculated effective detection volume is 3.2 m'- For this volume, a 
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concentration of 1 nM contends to an occupancy of 1.9. The difiusion time of a 
molecule of species a with a diffusion constant A. Anough this detection volume^ 
was calculated by simulating many paths through the detection volume: 
r"- (2.8x10- cm^)/z), . Since the persistence length of dsDNA is -150 bp or 500 
A [62], one can treat the DNA fiagments as rod-like polymers with length 
i«240A. and diameter 6«20A. with averaged translation^ difiusion constant 

7is the dynamic viscosity of the solution =i.i6 mPa s for aqueous soludon with 
5% v/v glycerol; [63]]. For r=25 -C, Z,=3.9xlO-' c^Vs, using this difiusion 
constant, the DNA fragu^cnts have a diffusion tin.e of .^=710^. in the calculated 
detection volume. We model the RNAP holoenzyme (dimensions 120x140x150 A. 
[57]) as a sphere with radius in the range 60-75 A, where 

D = k^/6x7}R. ^^32)) 
In this case, tiie e>q)ected difSision time is 0.9-1.1 ms. 

[0114] Fitting routine: The Levenberg-Marquardt nonlinear least square fitting 
procedure [59] is used to extract the parameters ftom the PAID histograms 
calculatmg the necessary partial derivatives numerically. To estimate the standard 
deviation of eirors, we use ten independent instances of the histogram to calculate an 
estimate of the statistical errors of each bin. The statistical errors are used as 
weights in the fitting routine [58]. If a bin is nonzero in fewer than ten instances it 
IS excluded ftom the fit. Without this restiiction, we found that only a few points 
dommate the value for (the^^ merit iunction is not ideal for sparsely populated 
bins in histograms). The ten independent instances were individually fitted (wifli a 
vanety of initial conditions, checking for convergence), extracting the means and 
standard erwrs of the model parameters. The standard enx>r of the mean with 10 
instances is ~3 where N^iois the sample size) times smaller than tire sample 
standard deviation for each parameter. The same procedmB is used for FIMDA, 
FCS. and FIDA. For each method, the model used was specialized to diffiisioii 
wathm a Gaussian detection volume (excluding triplet state fluctuations or other 
dynamics). FIMDA histograms were formed using the same bins as the PAID 
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histograms, and the spacing between counted inteivals was set to 10 ms for aU time 
bin widths. 

#115] For the experimental data sets, 300 s of data were split into 10 sections of 30 
s each; PAID histograms were calculated for each section, and the standani 
deviation of each bin was calculated to estimate the error of each bin; each section 
was fit separately, and the mean and error of the mean of each fitted parameter was 
calculated. Because the PAID model does not yet account for triplet state related 
fluctuations, the range of fitting for the experiments is restricted to time intervals 
greater than 10 jxs to minimize the effects on fitting. The procedure was modified 
for the low-occupancy dual-chamiel experimental data since there were too few 

bursts firom Dl^A^'"^^'^ and DMACys.iT ■ , 

^ ana UNA , resultmg m lai^ge errors for the 

extracted diffusion times. To obtain better statistics, die histogram for all 300 s was 
fitted, and subsequent bootstrap sampling was used to obtain error estimates for the 
extracted parameters [64]. To obtain one bootstrap sample, the data wer« spUt into 
10 30-s sections; 10 of these sections were ramlomly selected with replacement (die 
same section can be selected multiple times, or omitted). The PAID histograms for 
the selected sections were then averaged. We obtained 10 bootstrap samples in this 
mamier, and fit each of the resulting histograms; the standard deviation of the 10 
values extracted for each parameter is quoted as error bars for low-occupancy data. • 
I0116J Graphical representation of single-channel examples using PAH): The 
PAID histogram in accordance with the present invention is useful simultaneous 
measurements of diffusion time, brightness, and occupancy. We now describe how 
the parameters of a single diffusing species (occ^ancy, c, ; brightness per molecule 
in chamiel (.4=1), g„ • diffosion time, r," ; and background count rate, k,, ) aifect the 
PAID histogram for a single-channel experiment. Fig. 6A shows the PAID 
histogram for a 30 s simulation with o,=0.i. .."=,00^3, ,..=50kHz. and 
*b,=OkHz. At time interval ^=0. a photon fit,m one of the diffusing molecules is 
detected; fee position probability distribution for fee molecule that emitted the 
photon matches V^. At r=0. the autocorrelation amplihide from FCS is 
C«(r=0)=(c^+l)/ci [65]. Thismeansthattherear«,onaverage,c,+l molecules 
of species a =1 inside (since the molecules are assumed to be independent and 
there are nmnerous other molecules in a solution of a lar^e vohmie, the presence of 



wo 2004/011903 



PCT/US2003/023252 



-55- 



one molecuJe of a species does not affect tbe probability of a second molecule being 
present ). While the n>oIec«le that enutted the photon is present (time interval 
r rn, the average count rate in channels is {l+c.),„ . After a time interval 
r, T-. the molecule that emitted the photon has difiused out of v^, and the 
averagecountrateretumstothevalueforanarbitrarytime, There is 1 ridge of 
high correlation density at small . , (Fig. 6 A, y.hite and red contours) that decays 
wrth a time scale . (as with autocorxeladon in FCS). Since the count rate at 
this time interval is the peak of high correlation density follows a 

taiectory m the histogram of the fom. logn, .og[(l.,)...]..og.. wher. 
MO+^-kj the "vertical oflfeefofthe red ridge in Fig. 6A. After tire molecule 
diffi-es out of , the peak decays in height, indicating that the molecules present 
are unco^elated with the initial photon. The count rate decreases to tire average 
count rate. . and the correlation density peak follows a trajectory in histogram of 
thefomr logn log (c.^,.)+,og.. where log(c,^,.) is the vertical offiet of tire second, 
orange ridge in Fig. 6A (jipper-right comer). 

[0117] individual parameters influence specific features in the PAID histograms 
(Frgs. 6B-6E). When occupancy incr^ses 10-fold (Fig. 6B). fte histogram becomes 
narrower along tire monitor photon count axis (shown by opposing arrays),- tire 
. correlation amplitude decrees 10-fold, and tire long time count mte increases 10- 
fold (shown by arro^ on tire upper-right). When tire brightness increases 10-fold 
(Frg. 60, tire histogram shife up in tire log-log plot (as indicated by tire arrets) 
and becomes narrower along fl:e monitor photon axis. When tire diffirsion time 
mcreases 10-fold (Fig. 6D). tire red and white contours are extended 10-fold along 
tire tme interval axis as well as the monitor photon count axis (shown by arro^) 
When a constant background witir rute = 5 kHz is added (Fig. 6E), tire correlation 
density decreases (since many start photons are now uncorrelated background 
photons), tire background component is seen as an additional shallow slope (shown 
by tire arroW) of tire histogram to tire right of tire mdn correlation peak, and tire total 
count rate is doubled (seen at long time intervals r). In Fig. 6F. a second brighter 
species witir,,.=,..x4, c.=0.006. =100^s,and A,. =0.765 kHz is added; tirese 
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parameteis were chosen so that the coirelation curves conesponding to Figs. 6A and 
6F would overlap, danonstrating the limitations of FCS. 

10118] The advantages of PAID over FCS can be seen in Fig. 6G. which shows the 
autocorrelations for each of the simulations in Figs. 6A-6F. In FCS, the occupancy 
is extracted &om the correlation amplitude, which is C«(r=0)=l+l/ci for a single 
diffusing species with no background Qn Fig. 6A, c, =0.1, so C^(r=0)=11.0 for 
black curve). One FCS limitation is that both occupancy increases and 
background increases decrease the correlation ampUtude [cf red curve (ci'=10xc) 
and cyan curve (widi =5 kHz) in Kg. 6G]. In contrast, with PAID, they are 
distinguishable (Figs. 6B and 6E). Another FCS limitation is that changes in 
molecular brightness do not affect the correlation curve [cf black cun>e and dashed 
green cu,-ve {q'u = lOx?,,) in Fig. 6F], preventing the abihty to distinguish between 
components with different brightness (cf overlap of magenta line and black tine in 
Fig. 6G to difference between Fig. 6A and Fig. 6F). Again, the PAID histograms 
(Figs. 6A and 6C) can clearly distinguish such changes. 

[0119J Simulations-quantitative analysis of stoichiometry using single-channel 
PAID: To demonstrate the abihty of PAID to extract simultaneously occupancy, 
dififosion time, and brightness, and to compare it to FCS, FIDA/PCH, and FIMDA.' 
we generated two seri^ of single-channel (start = stop = monitor channel), single- 
component simulations in a Gaussian detection volume. The comparison of PAID 
with FIMDA is especially interesting, since both can simultaneously extract 
■occupancy, diffusion, and brightness. These simulations provide comparisons 
between the methods under representative conditions where aU methods have an 
appropriate model; we chose to use a Gaussian detection volume, ignoring triplet 
state induced flucniations. FIDA and FIMDA, PCH, and PAID account for non- 
gaussian detection volumes differently, and PAID does not yet account for triplet 
state induced fluctuations. The low-occupancy series consists of 10 simulations 
with c,=0.l, zf=100Ms, 9„=50kHz, and A„=0kH2. The intennediate- 
occupancy series is similar, but with c, =1.0. PAID perfonned well in botii series, 
with average error for each parameter <2% (Table 1); see Fig. 7 for a representative 
PAID fit for a low-occupancy simulation fix)m Table 1. For low-occupancy 
simulations, FIMDA extracts values with accuracy similar to PAID, but with a 
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significanfly worse fit (ar^~4; mainly for time delayxUffiision time [66]). For the 
intennediate-occupancy simulations, the fits are good (x'~l), and the extracted 
panunetere are close to the simulation values (except for a 5% downward bias in the 
difRjsion lime). The accuracy of FCS-extracted parameters was similar to PAID and 
FIMDA. The errors of the FEDA-extracted parameters were similar to the errors 
found using the other methods. The time bin width used for the simulations was 
chosen to be 1/5 of the shortest diflRision time (20 ^), since FIDA assumes the 
molecules are stationary during tire counting interval. Despite the short time bin, 
there is an upward bias in the occupancies, along with a downward bias in the 
brightness due to species difiusion during the time bin. The magnitude of these 
biases can be estimated using the corrections for the apparent occupancy 
'=m and brightness =qr^^ used in FIMDA [32]. is calculated for 

Gaussian detection volumes using y? = a>V/' =0.04 and t ^t/t" =0.2. We obtam 
=0.94 , which gives the expected values of =0.107 and 9^ =47 ktts, which 
match well with the values given in Table 1 . 
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Table 1: Parameters Extracted using PAID, FIMDA, FCS, and FIDA 
fits for single-channel, one-component simulations in a Gaussian 
detection volume 



Panuneten 


SimaJadon 


PAID 


FIMDA 


FCS 


FIDA 


Averaged Fits for 10 Simulations (30 s each): Low Occupancy 






1.2 ±0.1 


4.7 ±0.2 


1.3 ±0.1 


0.6 ±0.1 




0.0 


0.01 ±0.0] 


0.01 ±0.01 


N/A 


0.02 ±0.01 


C| (mol) 


OJ 


O.IOO±0.001 


0.100±0.00l 


0.099 ±0.001 


0.105 ±0.001 


(ps) 


100.0 


98±I 


99±] 


]00±1 


N/A 


^j, (kHz/hK)]) 


50.0 


48.9±0.3 


49.7 ±0.3 


N/A 


47.3 ±0.3 




Averaged Fits for 10 Simulations (10 


s each): Intennediate Occupancy 








1.0 ±0.1 


1.0 ±0.1 


1.0 ±0.1 


0.7±0.1 




0.0 


0.03 ±0.01 


0.1 3 ±0.03 


N/A 


0.3 ±0.1 


C| (mo]) 


1.0 


0.99 ±0.01 


0.99 ±0.01 


0.99 ±0.01 


1.05 ±0.01 




100.0 


99±1 


94±l 


100 ± I 


N/A 


(kHzAnol) 


50.0 


49J±0u! 


50.8 ±0.2 


N/A 


47.6 ±0.2 



Fixed values are in italics, fitted parameters are quoted with error estimates 
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(01201 We also perfonned low-occupancy (c=c,= 0.05) and intennediate- 
occupancy = 0.5 and ^ =c= 2.5) simulations to test PAID for detection of 

stoichiometiy in two-component experiments (Table 2). PAID perfonned well- low- 
occupancy simulations yielded low enois (lo/„.2-/o); intennediate-occupancy 
simulations (^=c,= 0.5 ) yielded somewhat higher enoxs (l%-8o/o). as weU as 
biases up to 10o/„ for the diffusion time and occupancy of the dimmer species For 
somewhat higher occupancy =^= 2.5). the quoted enors are still <lOo/.. but the 
biases reach up to 20o/o. The enors and biases of the FIMDA-extmcted parameters 
were m general similar to those extracted using PAID. However, for low-occupancy 
smiulations, the occupancy biases were larger (~20o/,). For intennediate-occupancy 
simulations (c, =c, = 0.5), occupancy biases were reduced, but the background count 
rate values were off by 25%. For the higher occupancy simulations (c, =c, = 2.5). 

^'^o^ were in the range of 10-300/0, with biases similar to PAID (up to 20O/O) ITxe' 
differences between PAID and FIMDA at low occupancy are primarily related to the 
models used, rather than flie histograms used [66]. FCS fits with fixed brightness 
drd not conver:ge for either high- or low-occupancy simulations; significantly 
drfiferent results with the same were found with the same data set This problem 
arises because FCS relies on difiusion time to detect the presence of subpopulations 
If both occupancy and brightness are fixed, difiusion thnes can often be extracted, 
but with poorer accuracy than with PAID or FIMDA; since FIDA/PCH can (at best) 
fix die occupancies and brightnesses to their conect values, this demonstrates clearly 
that simultaneous fitting of FCS and FIDA/PCH camiot match the perfonnance of 
PAID or FIMDA. The errors of the parameter, extracted fiom the two-species 
sunulatrons were larger for FIDA than for PAID or FIMDA. 
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Table 2: Parameters Extracted using PAID, FIMDA, PCS, and FIDA fits 
tor single-channel, two-component simulations in a Gaussian detection 
volume 



Parameters 


Simalation 


PAD) 


FIMDA 


FCS 


FIDA 




Averaged Fits for 10 Simulations (30 s each): Low Occupancy 








0.77 ±0.04 


1.1 ±0.3 


I.3±0.2 


0.5 ±0.1 




2M 


2.00 ±0.01 


].90±0.0] 


2.0 


1.99 ±0.03 


^ (inol) 


0.05 


0.048 ±0.001 


0.064±0.001 


0.05 


0.059 ±0.002 




100.0 


100±2 


117±4 


105± 7 


N/A 


^,1 (kHz/mol) 


mo 


50±1 


51 ±2 


50.0 


48 ±2 


C2 (mol) 


0.05 


0.052 ±0.001 


0.042 ±0.002 


0.05 


0.048 ±0.003 




150.0 


145±2 


146 ±3 


142±3 


N/A 


•^21 (kH2/mol) 


100.0 


98±1 


106±I 


100.0 


100±2 


Averaged Fits for 10 SimuJations (JO s each): Inteimcdiale Occupancy 






0.9 ±0.} 


0J5±0.02 


].4±0.2 


0.7 ±0.1 




2.0 


2.03 ± 0.02 


1.6 ±0.2 


2.0 


].6±0.3 


C, (mol) 


0.5 


0.45 ±0.03 


0.51 ±0.04 


0.5 


0.57 ±0.03 


r,°(ps) 


100.0 


88 ±6 


92±5 


120 ± 7 


N/A 


(kHz/knoI) 


50.0 


53±1 


52±2 


50.0 


46±4 


C2 (mol) 


0,5 


0J3±0.02 


0.50 ±0.04 


0.5 


0.52 ±0.04 


.^2° (MS) 


150.0 


153 ±1 


144±4 


141 ±3 


N/A 


^2] (kHz/mo!) 


100.0 


96±] 


101 ±2 


100.0 


97±2 




Averaged Fits for 10 Simulations (10 


s each): Intenxiediate Occupancy 








0.63 ±0.03 


0.5 ±0,1 


0.5 ±0.1 


0.5 ± 0.1 


^01 (kHz) 


2.0 


2,3 ±0.2 


4±2 


2.0 


0.4 ±0.4 


Cj (mol) 


2.5 


2.8 ±0.2 


2.7 ±0.3 


2.5 


3.0 ±0.3 


7 ° (MS) 


100.0 


1I0±I0 


100±30 


t27± 6 


N/A 


g,i (kHz^moO 


50.0 


55 ±3 


49±2 


50.0 


50±4 


C2 (mol) 


2.5 


2.1 ±0J 


2.2 ±0.2 


2.5 


2.2 ±0.3 


^2 (MS) 


150.0 


150±8 


143 ±5 


136± 3 


N/A 


^21 (kHz/moI) 


100.0 


101 ±2 


108± 3 


100.0 


I02±5 
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[0121] Experiments-quantitative analysis using single-channel PAID: 
Measurements were perfonned on DNA fiagments to test the ability of PAID to 
detect subpopulations in solution based on the properties available in a single 
channel. Three samples were tested in the low-occupancy regime: 0 1 nM 

DNA^^'^B.'T only; 0.1 nM DNA^>'^b..ta:v3b.«b ^ ^ ^ 

Dl^A^y^^-^^ and 0.05 nM DNA^3B.,t«:.3b.65b (tables 3-4). Similarly, three samples 
were tested in the intennediate-occupancy regime: 1 nM DNA^^-'^ only 1 nM 
only; and a mixture of 0.5 nM DNA^^^^-'^ and 0 5 nM 
uiNA (Tables 3-4). For the fits of the one-species samples listed in 

Table 3, all parameters were fitted; the model fits well (x' in the range of 0.8-1.5). 
The important feature of the analysis is that the dnA'^^^^-''^'^''^^'*'^ is ~2.1 times as 
bright as the DNA«=y».'T (Table 3). demonstrating die ability of PAID to evaluate 
stoichiometry; similar stoichiomeHy results are obtained using PCS in combination 
with the mean count rate. Based on the calculated detection vohmie, the 
occupancies of 1 nM and 0.1 nM samples are expected to be, respectively. 1.9 and 
0.19. The occupancies extracted are wifliin 30% of these values (enora attributed to 
pipetting errors and losses on surfaces). The measured difiusion times (620-760 ys) 
match well to the calculated diffusion time (710 ps). The measured background 
rates for three of the samples (0.2-0.4 kHz) are consistent with buffer-only 
measurements (0.24 kHz); the exception is die background in the intennediate- 
occupancy sample of DNA^^«.'T«^».«B; however, the background makes up <lo/o 
of the signal. 
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Table 3: Parameters Extracted using single-component PAID fits for 
smgle-channel, single-species experiments 



Parameters 






Averaged Fits foi 


10 Measurements (30 


s each): Low Occupancy 




0.8 ±0.] 


1J±0.1 


*bl <'^) 


039 ±0.03 


0.37 ±0.04 


<\ (mol) 


0.1 9 ±0.0] 


0.14 ±0.0] 


r,°(Ms) 


640 ±20 


760 ±30 


qu (kHz/inol) 


9,9 ±0J 


21.5 ±0.7 


Averagea i-its tor 10 Measurements (30 s each): Interaiediate 




1.5±0.1 


I^±0.1 


K\ o^Hz) 


0,23 ±0.04 


2.2 ±0.3 


Cj (mol) 


2.21 ±0.02 


2.32 ±0.04 


7° Ok) 


620± 10 


650± 10 


^„ (kHz/moI) 


8.7 ±0.1 


1 7.9 ±0.2 



Buffer-only measurements - k„ = 0.24 ± O.OI kHz 
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Table 4: Parameters Extracted using two-component PAID fite 
smgle-channel, one- and two-species experiments 



C, (raol) 

9i, (]cH2/jnol) 
<moI) 

^-2° (US) 

^21 (kHz/moI) 



Restricted Fixed ratio fit Restricted fit DNA«wn«*«M5B Dj^j^cmiMwi^i 

■ , Kxedrattj fit Restricted fif 

Averaged Fib for 10 M easurements (30 s each): Low Occupancy 

2.3±0J 1.1 ±0.1 

OM 0.40±0.Q2 

0.01 ±0.01 0.10 ±0.02 
640 760 ±20 

^9 9.6 ±0.4 

O.I 6 ±0.01 0.08 ±0.01 

^0 760 ±20 

^^'^ 19.3±0.8 
Averaged Hts for 10 Measunanents (30 s each): Intoinediate Occupancy 





0.8±0.1 


1.9±03 


1.2±0.1 




036 ±0.02 


0.24 


035 ±0.03 


Cj (mol) 


0.1 9 ±0.01 


0.210 ±0.002 


0.01 ±0.01 


(US) 


680 ±10 


640 


760 ±20 


9, J (kHz/niol) 


9.7 ±0.2 


9.9 


1 0.8 ±03 


^2 (mol) 


0,01 ±0.01 


0.001 ±0.001 


0.14 ±0.01 




680 ± 10 


640 


760 ±20 


^21 (kHz/moI) 


4,8 ±0.1 


19.9 


21.6±0J 



1.9 ±03 
0^4 
0.1 30 ±0.004 
640 

9.9 

* 0.068 ±0.003 
640 
J9.9 



13±0.1 


1.6±0.1 


1.0±0.I 


13 ±0.1 


2,5 ±0.1 


2.8 ±0.1 


0.7 ±0.1 


0.24 


0J±0.1 


0Ji4 


0.9 ±0.1 


0.24 


2.04 ±0.04 


2.19 ±0.01 


0J±0.1 


0.2 ±0.1 


1.4 ±0.1 


1.99 ±0.04 


640±]0 


620 


630 ±10 


620 


680 ± 10 


620 


9.0 ±0.2 


8.7 


92 ±0.1 


8.7 


7,4±0.1 


8.7 


0.10 ±0.03 


0.005 ±0.002 


2.1 ±0,1, 


2.41 ±0.04 


1,0 ±0.1 


0.49 ±0.03 


640+10 


620 


630 ±10 


620 


680±10 


620 


4.5 ±0.1 


17.4 


18.5 ±03 


17.4 


14.9 ±03 


17.4 



difiusion times of flie two components are linked. 

* The brightness values are 1 and 2 times the values extracted from the DNA'^'^-''^ 
samples, and the diffusion times are fixed. The background rates were extracted 
from separate experiments. 
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, I0122J PAID examples were performed using Cy3 (rather than Cy3B); in this case, 
the ratio in brightness between DNA'=y^'Tycy3.fisB jy^^C3,n ^ 
shown). Even at low excitation intensities (<10 kW/cm^ ), the factor of -1.5 
persisted, excluding triplet state saturation and photobleaching from tiie possible 
sources of the discrepancy. Using FCS. we observed a fluctuation that became 
faster with increasing excitation intensity without a corresponding change in the 
fluctuation ampUtude (unlike triplet fluctuations). A similar effect was observed 
previously for Cy5, and was identified as photo-induced isomerization [67]. Since 
isomerizalion of two distinct fluorophores is uncoirelated, the amplitude of the 
• fluctuation is reduced by a fector of 2 in the DNA^^-'^'^^.ese^ increasing the 
apparent concentration of the double-labeled species as compared to the single- 
labeled DNA^^'T. resulting in only a fector of -1.5 difference in brightness. Using 
Cy3B. a confonnationally-constrained analog of Cy3 which prevents isomerization, 
the fluctuation is absent. This emphasizes the importance of the choice of 
fluorophore when attempting to use brightness as a measure of stoiohiometry. Dual- 
channel methods (discussed below) are less sensitive to such effects. 
[0123] To determine if two species are present in a sample, it is necessary to 
perform two-component fits as weU as one-component fits. Table 4 hsts the results 
for two-component fits of the one- and two-species samples. For the "fixed ratio 
fit", the brightness of one species is fixed to be twice the brightness of the other, the 
total brightness can vary and the diffusion times of the two components are liriked 
(set equal to each other). For the "restricted fit", only the occupancies are aUowed to 
vary; the brightness values for the two components are set to be 1 and 2 times the 
brightness fiom the single-labeled species in Table 3. the diffusion times are taken 
from the same species, and the background is taken from separate, buffer-only 
measurements (not shown). For the single-species samples, only one component is 
fitted with a significant occupancy for most of the fits (ratio of occupancies >10:1); 
the °°ly^^^^'^eption is the "fixed ratio fif of the intemiediate occupancy sample of 
jy^j>^Cy3B.n/Cy3BfiSB ^^^^^ ^.j. j^^^^^^^^ restricted fit has a ratio of 12:1). For 
the mixture samples, significant occupancies are fitted for both components for all of 
the fits. By dividing the single-species occupancies (extracted in Table 3) by 2. we 
obtain the expected occupancies for the mixture; the expected occupanciel of 
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DNA<^3B.,T DNA<^3B..T«>,3B^ ^ respectively 0.10 and 0.07 for the low- 
occupancy samples and 1.1 and 1.2 for the intermediate occupancy samples. The 
occupancies extracted from the low-occupancy mixture match the expected 
occupancies wifliin 10% (except the lower brightness species is fit with a 30% 
higher occupancy with the restricted fit). The occupancies extracted 6om the 
intermediate occupancy mixture are biased toward the lower brightness species 
(identified wifli DNA^^^-'^; 30% and 80% higher in the fixed ratio fit and restricted 
fit, respectively), decreasing the amount detected in the higher brightness species 
(identified with DNA^^B-' W5b. ^qo/, i^^,, 50% lower in the fixed ratio fit 
and the restricted fit, respectively). These biases are most likely due to inadequate 
modeling of the detection vohune; spurious components or biases are generally the 
result of imperfect overlap of model and data [68]. Improvements are expected 
upon measurement of the confocal detection volume and modeling of additional 
photophysical properties of fluorophores. Nevertheless, significant occupancies for 
two components were extracted in the mixtures, but not in the single-species 
samples; this demonstrates the capabilities of PAID to detect heterogeneity in single- 
channel experiments. For flie fixed ratio fits, the diffusion times extracted (630-760 
MS ) match exactly the calculated value (710 ixs). The brightness -values are 
consistent with the values extracted fi-om the single-component fits (aU are within 



[0124] Graphical representation of dual-channel experiments using PAID: 
When appUed to dual-channel examples the ability of PAID to extract brightness in 
multiple channels for each species becomes extremely inqjortant. Accordingly, fee 
use of PAID for dual-channel systems is preferred. We use the example of a simple 
binding assay to illustrate what to expect when using PAID in dual-channel 
experiments. As illustrated in Fig. IB, in a dual-color fluorescence bmding assay, 
one molecule A^ is labeled with a fluorophore of one color (for example "yellow", 
denoted:,;), while the second molecule B' is labeled with a fluorophore of a second 
color (for example "red", denoted r). Each fluorophore is excited by a distinct laser 
wavelength. Upon binding, diere are three species present in solution: fi-ee A-*- , free 
B^ and complexes A^^ The "yellow" and "red" fluorophores y and r have 
corresponding "yellow" and "red" detection channels, denoted, respectively, Y and 
R. Because of the vibronic tail of organic fluorophores toward the red end of their 
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emission spectra, there is a smaU contribution &om y into tbe channel R; the 
contribution firom r into channel 7is typically negligible. 

[0125] For each two-channel data set (firom experiment or simulation), we 
perfonned a global fit of the series of all possible dual-channel PAID histograms, 
extracting brightness (in both channels), difiiision time, and occupancy for the 
species A^^ , A" and B^. We write the channel assignments as a three letter code. 
STM, specifying the start, stop, and monitor channels (eg. the RYR PAID histogram 
uses S=R, T=Y, and M=7?). For dual-channel experiments, there are 8 unique 
channel assignments for the PAID histogram: RRR, RRY, RYR, RYY, YRR, YRY, 
ITR. and 3T7. The assignments of the start and stop channels (firet two letters) 
select species of interest, which emit photons in both the start and stop channels. 
The choice of monitor channel detemiines the fluorophore which wiU be analyzed in 
teims of brightness using the PAID histogram. Figs. 8A, 8D. 8G. and 8J show 
histograms for a dual-channel experiment if background, A^* . B"- , and A'B' are 
present. Figs. 83, 8E. 8H. and 8K show the histograms if A^' is absent. The 
cartoons of the free molecules and complexes indicate which species contribute to 
any large correlation density peak. Figs. 8C, 8F, 81, and 8L compare vertical sUces 
at time interval r = 1 ms for the histograms in the presence and absence of A'B' . 
10126] Taken together, the RYR, RYY, YRR, and YRY histograms detennine the 
occupancy, difiusion time, and brightness in the 7and i? channels of A^' (these are 
direct extensions of the cross-correlation functions used in FCCS). The histograms 
with S=R and T=Y emphasize time regions when A^' is present, smce A^' emits 
photons in both i? and 7 whereas A' and B' emit only in one channel. SpecificaUy. 
the RYR histogram monitors the brightness of r (Fig. 8A: RYR with A^' present; 
Fig. 8B: RYR with no A^O- The correlation density peak in Fig. 8A corresponding 
to A^' (shown by arrow) is absent in Fig. 8B (flie small correlation peak in Fig. 8B 
corresponds to the contribution of j; into R). The RYY histogram (not shown) 
monitors the brightness ofy, and is similar to RYR. The YRR and »?r histograms 
(not shown) are also similar to the RYR histogram in this case. 
[0127] The histograms with S=R and T=R emphasize time regions where B' and 
A^' are present, since both emit photons in R. Specifically, for the RRR histogram, 
the correlation peaks fiom B'and A^'overlap since both have a similar brightness 
in R (Figs. 8D-F). There is only one correlation peak visible, where B' and A^B' 



wo 2004/011903 



PCT/US2003/023252 



-67- 



both contribute (Fig. 8D, arrows), so &e RRR histognun cannot distinguish weU 
between and A^. However, for the iWF histogram, the coirelation peaks &om 
B' and A'B' aie well separated. A'B' emits in 7 (resulting in a correlation peak 
with high monitor photon count; left arrow), whereas BT does not emit in 7 
(resulting in a correlation peak with low monitor photon counts; right arrow) (Figs. 
8G and 81). When is absent, the coiresponding coirelation peak is noticeably 
absent, leaving only the coirelation peak resulting from B' (Figs. 8H and 81). The 
ITT and YYR PAID histograms emphasize A^-and A^'in a similar manner (Fig 8J- 
L). 

[0128J Higher sensitivity could be achieved if two monitor channels were used for 
the two-channel experiments, allowing the use of ratiometric infoimation not 
available with a single monitor channel. Observables in single-molecule difiusion 
studies that depend on the ratio of two channels are more sensitive than brightness. 
A PAID histogram with two monitor channels has three axes: the time interval axis 
r . the number of photons in the red channel n,, and the number of photons counted 
in the yellow channel n, . Fig. 9A shows three shoes of this histogram at different 
time intervals when A'. B', and A'B' are present, and Fig. 9B shows the same 
three slices when A'B' is absent. The histograms shown are for the PAID 
histogram with the red channel assigned as the start channel (S = R), the yellow 
channel assigned as the stop channel (r = 7). and a monitor channel for each 
chamiel (M, = R and M, = 7). A cartoon of a particular type of molecule is placed 
in the r = 1ms slice at the approximate location where it contributes. The 
histograms in Figs. 8A and 8B are the collapse of Figs. 9A and 9B, respectively, 
summing along the yellow monitor channel axis, while keeping the time interval and 
red monitor channels axes. For an isolated burst, the number of photons counted 
depends on the difiusion path taken dirough the detection volume, whereas the ratio 
between two channels does not For example, the coirelation density peak in Fig. 
9A at T = 1 ms corresponding the complex is wider in the direction of the diagonal 

=n, than in the perpendicular direction. The one monitor channel PAID ' 
histogram is a collapse of the two monitor channel histogram onto the 7 or /J axis. 
This collapse smears the central peak; the subpopulations of complex and free 
molecules are not cleariy separated in the one monitor channel case, but are 
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separated with two monitor channels (compare the r =1 ms sUce in Fig. 9A with the 
coiresponding single monitor channel histogram in Fig. 8A.) In view of these 
benefits, the fitting model can be extended to account for two monitor chamiels. if 
desired. TheremainderoftheexampleswiHuseonlyasinglemonitorchamiel. 
[0129] Simulations-quantitative analysis of binding using dual^Aannel PAID: 
PAID and FCCS in dual-channel applications with two spectrally-separable 
fluorophores were investigated using multiple species simulations. The diffosion 
time of A'-B' was set to be 33% higher than A" and B- the biightiiess of each 
fluorophore was unchanged. As in Ae single-channel case, low- and intermediate- 
occupancy reghnes were studied (Tables 5 and 6). The fluorophore y was set to have 
a total brightness of 50 kHz. 90% in chamiel 7 and IQo/o in chamiel R. The 
fluorophore rwassetto haveatotalbrightaess of SOkHz. 100% in channel R, The 
background was set to =k^g=2 kHz . 
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Table 5: Parameters Extracted using PAID and FCCS fits for two- 
channel, three-component simulations in a Gaussian detection volume at 
low occupancy 



Parameten 


Simiilatkrii 


PAID- 
Unrestrf cted Fit 


PA10- 
Rcstricted Fit 


FCCS 


Averaged Fits for 1 0 SimuJations (30 s each): Low Occupancy 


%^ 




0.72 ±0.02 


0.79 ±0.02 


0.95 ±0.04 


Kr (kHz) 


2.0 


2.00 ±0.01 


2,0 


2.0 


Kr CkHz) 


2,0 


2.00 ±0.01 


2.0 


2.0 


Cj (mol) 


0.05 


0.049 ±0.001 


0.049 ±0.001 


0.050 ±0.001 




300,0 


287 ±7 


300.0 


31I±4 


q^g (kHzAnoI) 


50.0 


50.3 ±0.4 


50.0 


50.0 


qiy (kHz/mol) 


0.0 


0.05 ±0.01 


0.0 


0.0 


C2 (mol) 


0,05 


0.049 ±0.001 


0.050 ±0.001 


0.050 ±0.001 


(Its) 


300.0 


294 ±7 


3oao 


306 ±4 


q^n (kHz/moI) 


5.0 


4.8 ±0.1 


5.0 


5.0 


qjY (kHz/mo!) 


45.0 


44.8 + 0.3 


45.0 


45.0 


Cj (mol) 


0,05 


0.050 ±0.0Q) 


0.050 ±0.001 


0.051 ±0.001 




400.0 


394±6 


389 ±7 


404 ±8 


^3jj (kHz/mol) 


55.0 


54.9 ±0.3 


55.2 ±0.1 


55.0 


q^Y (kHz/bo!) • 


45.0 


44.7 ±02 


44.8 ±0.1 


45.0 



Values that are fixed are shown in italics, with no enors listed. 
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Table 6: Parameters Extracted using PAID and FCCS fits for two- 
channel, tbree-component simulations in a Gaussian detection volume at 
intermediate occupancy 



Parameters 


Simulation PAII>- 

Unrestricted Fit 


PAID. 
Restricted Fit 


FCCS 


Averaged Fits fc 


r to Simulations (10 s each): bitemediate Occupancy 








0.81 ±0.02 


0.90 ±0.04 


0.76±0.09 


Kr (kHz) 


2.0 


2.02 ±0.02 


2.0 


2.0 


K7 (kHz) 


Z0 


2.00 ±0.02 


2.0 


2.0 


C, (moO 


ftj 


0.50 ±0.01 


0.49±0.01 


0.050 ±0.01 


7° (MS) 


300.0 


320±9 


300.0 


3I3±3 


OcHz/moI) 


50,0 


49.1 ±0.2 


50.0 


50.0 


q^Y (kHz/moI) 


0.0 


0.06 ±0.02 


0.0 


0.0 


^2 (mol) 


0.5 


0.49 ±0.01 


0.49 ±0.01 


0.50 ±0.01 


(MS) 


300.0 


306±11 


300.0 


303 + 7 


^jJt, (kHzf^mol) 


5.0 


4.9 ±0.2 


5.0 


5.0 


OcHz/moI) 


45.0 


44.5 ±0.4 


45.0 


45.0 


C3 (mol) 


0.5 


0.50 ±0.01 


0.49 ±0.01 


0.49 ±0.01 


^3^^ (MS) 


400.0 


395 ±4 


390±6 


407 ±5 


q^ft (IcHz/moO 


55.0 


54.3 ±0.4 


55.0 ±0.2 


55.0 


^3,. (kHzAmoI) 


45.0 


44.5 ±0.5 


45.1 ±0.2 


45.0 
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[0130] For each simuJation, three fits were perfoimed. The l** fit (Tables 5 and 6, 
3"» column) uses PAID with all parameters unrestricted The 2'"' fit (Tables 5 and 6, 
4* cohmin) also uses PAID, with all parameters of the ftee components except for 
the occupancies being fixed. For the l" and 2°^ fits, all combinations of the dual- 
channel PAID histogram are fitted simultaneously; since each histogram emphasizes 
different species, parameters for all species can be extracted. The 3"* fit (Tables 5 
and 6. 5* column) uses FCCS to simultaneously fit the autocorrelations of the red 
and yellow channels and the two cross-correlations, with all brightness and 
background values fixed. 

[0131] The um-estricted fit extracted reliable values, within 1-10% for all parametei^ 
in both sets of simulations. Unexpectedly, fixing the brightness of the fiee 
components and the background resulted only in a modest improvement for the 
brightness values for the complex and no improvement for the difiusion time or for 
the occupancies. The parameters extracted usmg FCCS had similar statistical 
accuracy (-]•/, range) to those found using PAID. The extracted diffusion times are 
somewhat better with FCCS (the eiror bar^ were -20% smaller), mainly due to the 
smaUer numbers of fitted parameters (brightness and background were fixed to their 
simulatiori values; for experiments, it would be necessary to measure the brightness 
and badcground vahies using a method different than FCS.) 

[0132] Experiments-quantitative analysis using dual-channel PAID: We 
performed measurements on fluorescently labeled DNA fiagments to test the abihty 
of dual-channel PAID to detect multiple species in sohition and measure their 
properties, as is necessary for analysis of interactions. The fluoixjphores used in 
these experiments were Cy3 (as the «Vellow"fluorophorej;) and Cy5 (as the 
"red" fluorophore r). Based on the occupancies measured at higher concentration 
and the dilutions used, the occupancies for the intennediate-occupancy samples were 
expected to be 0.85 + 0.09 for DNA^-'^, 1.15 ± 0.05 for DNA^^^^-^^, and 0.84 ± 
0.03 for DNA^5,iTm65B ^j^^ fragments were prepared as free components and in 
the following mixtures: DNAC^^-'^/DNA'=>'^.«« (to simulate non-interacting species); 

and DNA*^''""^^'"^ /DNA*^^'''r/nMACy3.<5B 

.^i^j^ /ijiNA /DNA ' (to simulate mteractmg species). The 

same samples were also prepared for 10-fold lower occupancies (low-occupancy 

samples). The value for DNAC>W.65b was found by analyzing fte Cy3 

fluorescence, so the occupancy quoted includes DNA*'*'''''^'^^^ and 
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DNA^^'''"'*^''"^^'"® (with non-fluorescent Cy5 - as seen in previous studies [12, 
69]). 

[0133] Fig. 10 shows how to visually detect &e presence of the double-labeled 
DNA*^^^*''^^^^'^^^ using PAID histograms; two low-occupancy data sets are shown: 
DNA^^^'^^^^^.^^« /DNA^^'»^/DNA^^-«« (Figs. lOA, lOD, lOG, and lOJ) and 
DNA^y^'^^/DNA^^-^^^ (Figs. lOB, lOE, lOH, and lOK). Figs. IOC, lOF, 101, and 
lOL compare vertical slices at time interval r = lmsfor the two data sets. By 
comparing the two sets of histograms, we can identify jy^/iOy5,\T/cy2fi5B 
histogram (Figs. lOA-C) emphasizes dNA^^'^*^^^^'"^; without DNA^^^'^"^^^^'^^^, 
the peak shown by the arrow in Fig. lOA disappears. The RRR histogram (Figs. 
lOD-F) emphasizes DNA^'^^ and dNA^''^^^«^^^ but does not distinguish weU 
between diese species since they both emit equaUy in R ; absence of 
DNA^^'^*^'^^'^^^ decreases amplitude of the correlation density peak (cf. Figs. lOD 
and lOE). The RRY histogram (Figs. lOG-I) also emphasizes DNA^^'^-''^ and 
DNA^^^*^"^^^^*^^^, but is able to distinguish between them since DNA^'^*^ does not 
emit in 7 whereas BNA^^^'^^^^^-^^^ does. Finally, the YYR histogram (Figs. lOJ-L) 
emphasizes DNA^^^'^^« and DNA^^'^^^^^'"^ separating them by their emission in 
R. In this case, because the occupancy of DNA^^^-^^® is higher daan 
j)^ACy5,iT/Cy3.65B^ the peak corresponding to DNA^^^^-^^^^'^^^ is less pronounced. 
[0134] Using PAID, occupancy, diffusion time, and brightness was extracted jfrom 
dual-channel experiments performed on single-species samples (Tables 7 and 8), 
then from dual-channel experiments performed on mixtures of those species (Tables 
9 and 10). Consistency between tiie values extracted from fits of data from single 
species samples and values extracted from multiple-species fits of mixtures 
demonstrates the ability of PAID to analyze mixtures of species. 
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Table 7: Parameters extracted using PAID fits for two-channel, single- 
species, low occupancy experiments 



Parameters 








Fit for 1 Measurement* (300 s): Low occupancy 






5.5 ±0.1 


1.9 ±0.3 


(kHz) 


0.82 ±0.01 


0.76 ±0.01 


0.68 ±0.01 


Aor (kHz) 


1.30±0.0I 


1.14 ±0.03 


0.95 ±0.01 


DNA Fragment 2 








Cj (mol) 


0.022 ±0.001 


0.073 ±0.003 


0.023 ±0.003 




430±10 


570 ±10 


670 ±60 


qijf (kHzAnoI) 


12.0 ±0.1 


0.94 ±0.03 


0.8 ±0.1 


q^Y (kHz^ol) 


0.02 ±0.01 


10.0 ±0.3 


9.2 ±0.8 


DAA Fragmeail 


JVone 


None 




C2 (moO 






0.012 ±0.001 








650 ±10 


^2i? (kHz/moI) 






6.8 ±0.3 


^27 (kHz^O 






6.7 ±0.3 



*A11 300 s of low occupancy data fit at once to improve statistics. 
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Table 8: Parameters extracted using PAID fits for two-channel, single- 
species, intermediate occupancy experiments 



Parameters DNA*^**'" ©jsj^c^wa j^^^cgsjmctxssB 



Averaged Fits for 10 Measucements (30 s): faitenoediate occupancy 





4.7 ±0^ 


2.2 ±0.1 


1.4±0.1 


(kHz) 


130 ±0.01 


0,99 + 0.01 


0.92 ±0.01 


k^y (kHz) 


1J3±0.01 


1.44 + 0.01 


1.29±0.01 


DNAFragmetttl 








<\ (moO 


0.68 J: 0.01 


1.08 ±0.01 


0J2±0.01 


(MS) 


390 ±10 


550 ±20 


700±10 


q^n (kHzAnoO 


9.4 ± 0.2 


0.82 ±0.01 


0.70 ±0.03 


q^Y (kHz/raol) 


0.02 ±0.01 


9.0 ±0.1 


10.7 ±0.2 


DNA Fragment 2 


None 


None 




Cj (mol) 






0.24 ±0.01 








580 ±10 


q^R (kHz/mo!) 






7.6 ±0.1 


q^f OcHz/moI) 






7.6±0.2 



) 
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Table 9: Parameteis extracted using PAID , fits for two-channel, 
multiple-species, low occupancy experiments 



Parameters 



Unrestricted 
2 component 



Restricted 
3 component 



Unrestricted 
3 component 



DNA^"****», 

Restricted 
_ 3 component 



Fit for I Meastiiement (300 s): Low ocaqiaacy 





2.8 ±0^ 


3.5 ± O.J 


23 ±0.4 


3.9±0.1 




0.87 ±0.01 


0.88 ±0.01 


0.80 ±0.01 


0.78 ±0.01 




1.04 ±0.01 


1.04 ±0.01 


1.12 ±0.01 


1.08 ±0.01 


DNA Fragment! 










C, (mol) 


0.018 ±0.001 


0.01 6 ±0.001 


O.OII±0.001 


0.01 6 ±0.001 




430 ±10 


430 


360 ±10 


420 


(kHz/mol) 


11.0 ±0.2 


12.0 


15.4 ±0.4 


J2.5 


(kHz/mol) 


0.0 1 ±0.02 


0.02 


0.02 ±0.03 


0.0 


DNA Fragment 2 










Cj (mol) 


0.068 ±0.001 


0.067 ±0.001 


0.077 ±0.003 


0.088 ±0X01 


^■2^(MS) 


530±10 


570 


570 ±10 


570 


(kHz/mol) 


0.95 ±0.02 


0.94 


0.89 ±0.02 


0,94 


^27 (kHz/mol) 


10.2 ±0.1 


JQ.0 


n.2±0.3 


J0.0 


Z>A'/4 Fragments 


None 








Cj (mol) 




0.0001 ±0.0001 


0.015 ±0.002 


0.013 ±0.001 






650 


680 ±40 


650 


^3;^ (kHz/moO 




6,8 


7.3 ±0.7 


6,8 


q^y (kHz/biol) 




6.7 


6.8 ±0.3 


6.7 



' Includes contributions from DNA^^"^*^' 
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Table 10: Parameters extracted using PAID fits for two-channel, 
multiple-species, intermediate occupancy experiments 



Parameters 



Unrestricted 
2 component 



Restricted 
3 component 



DNA***»», DNA'^wsi., 

DNA***", DNA*^, 

Unrestricted Restricted 

3 component 3 component 



Averaged Fits fig 10 Measunanents (30 s cadi): Intennediate occupancy 





K8±0.I 


3.0 ±0.2 


3.0 ±0.1 


4.8 ±0.2 




1.4±0.1 




1.66 ±0.01 


0.7J 


Ky O^Hz) 


1.44 ±0.01 


US 


I3±0.1 


1.15 


DNA Fragment 1 




DNA^ 


DNA^ 




Cj (mol) 


0.43 ±0.01 


0.57 + 0.01 


0J7±0.01 


0.74 ±0.01 


(jis) 


3SO±]0 


389 


360±10 


3S9 


(kHz/mol) 


1].0±0.1 


9.4 


14.5 ±0-4 


9.4 


(kHz/mol) 


0.20 ±0.03 


0.02 


0.3 ±0.1 


0.02 


DNA Fragment 2 


DNA"^*" 




DNA<^* 




(mol) 


0.86 ±0.01 


0.96±0.01 


1.0±0.l 


1.56 ±0.01 


(MS) 


560±10 


554 


500± 10 


554 


(kHz/moI) 


0.75 ±0.02 


0.82 


0.5 ±0.1 


0.82 


(kHzAnol) 


9.8 ±0.1 


9.0 


11.5±0-3 


9.0 


iW/l Fragments 










£3 (mol) 




0.01 ±0.01 


0.58 ±0.05 


0.22 ±0.01 


^■3° (}iS) 




575 


450 ±10 


575 


^3;j (kHz/mol) 




7.6 


5.3 ±0.4 


7.6 


q^Y (kHz^ol) 




7,6 


6.4 ±0.1 


7,6 



* Includes contributions from DNA'^^'**^"'^*^' ^'^^^^'^^^ 
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I0135J Figure 11 shows a representative fit of the RYY histogram (emphasizing 
D|sjy^Cy5.!T/cv3,65Bj low-occupancy mixture sample (Table 9); this is one of the 

eight histograms fitted simultaneously over time intervals between 10 fis and 1' s. 
The peak toward higher photon counts in the vertical slice is caused by 
D>jA^5.jT/cy3,65B generally within the error bars, and foUows the data to 

within 10%. 

[0136] As shown in Tables 7-10, we obtained fits with x"- in the range 10>;jr= >!■ 
(ideaUy 1 ). The quality of the fits is sufficient to extract values with generaUy 
exceUent accuracy and consistency, but the deviations indicated by the higher x^ 
give rise to biases in certain values (see below). The deviations between model and 
fit can be attributed to two general sources. First, many photophysical properties of 
the fluorophores are not yet modeled in PAID, such as triplet state fluctuations and 
fluorescence saturation [15], photobleaching [16, 70], and photo-induced 
isomerization [67], Although time intervals less than 10 jis are excluded from the fit 
to minimize photophysical effects, these properties stiU affect fitted values. Second, 
deviations of the experimental detection volume fi-om Ae modeled volume may also 
play a role in the increased x^ [68]. The model used for the fits in Tables 7-10 was 
based on a non-Gaussian detection volume; using a simpler Gaussian detection 
volume, significantly higher;ir^ values are found. For example, fitting tiie 
m^pcy^m ^^^^^^ ^ ;ir'=32.5 for the lowH)ccupancy 

and intermediate-occupancy data sets, respectively; the non-Gaussian volume 
yielded^2=5.5 and 2.2 respectively. This highlights the critical role of 
detection volume for the model, and indicates that further improvement may be 
possible. 

[0137] Tables 7 and 8 show fitted values for the samples containing only one 
species of labeled DNA, where the background rates, occupancies, diffusion times, 
and brightness in each channel were fitted parameters. To fit the data for 
i^iNA , two components were necessary, one with Cy3-only and another 

with CyS and Cy5. This is due to non-fluorescent Cy5; single-species fits and 
inspection of time traces both indicate the presence of a species emitting only in Cy3 
channel [12, 69]. When parameters corresponding to DNA^^*'^ were added to 
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check if there was a species with non-emitting Cy3 (restricted to fit emission only in 
R), there was only a smaU occupancy fitted (-0.001, not shown). For the single- 
species, low-occupancy data (Table 7). the extracted occupancies were 35%-70% 
lower than expected from higher-concentration FCS experiments, whereas for the 
single-species, intennediate-occupancy data (Table 8), they are 4%-20% lower than 
expected (note: for DNA'^^-""'^^-"^ the occupancies for DNA^''^'^^-*'^ and 
jjjsj^Cys&rtiiTmfisB added). The difference between extracted and estimated 
occupancies can be attributed to loss of DNA on surfaces during handling. The 
difiiision times extracted for UNA'^''"'^'-^ and DNA*^^-*^ are similar to the 
theoretical ones (550-700 jis vs. 710 ps); diffusion times for DNA'^^-''^ are shorter 
(-400 ps), largely due to the effects of photo-induced isomerization [67], although 
photobleaching of Cy5 within the detection vohnne and triplet state fluctuations may 
play a role. The brightness vahies extracted are similar to what was found earlier in 
the single-channel experiments. 

[0138] We performed two fits for mixtures of DNA fragments simulatmg non- 
interacting and interacting species (Tables 9-10). The 1^' fit assumes the correct 
number of species, but allows all parameters to finely vary; the 2"' fit uses the 
single-species parameters akeady extracted to restrict the parameters for the free 
components, except for occupancy. These fits show that a sample with two species 
can be distinguished from a sample with three species, as is necessary for measuring 
interactions. In Tables 9 and 10, PAID performed extremely well in all the 
restricted fits and m fliree of the four unrestricted fits. Although the three- 
component, unrestricted fit at intermediate occupancy identified three species, it 
produced values less consistent with values obtained widi oflrer fits. 
[0139] The restricted fits m Table 10 fix die background rates to values obtained 
from buffer-only measurements. The fits in Table 9 aUow the background rates to 
freely vary since, under low-occupancy conditions, contributions from the 
background photons to the PAID histograms are well-separated from contributions 
from the DNA species (when the background values are fixed to those used in the 
restricted fits in Table 10, small variations in background between measurements 
cause to be as high as 40). The occupancies extracted from the mixtures using 
restricted fits were consistent with the occupancies extracted from the single-species 
samples (within 15%. except for DNA^^^''^ at low occupancy, which is 27% lower ). 
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[0140] For the unrestricted low-occupancy data fits in Table 9. the occupancies 
extracted for DNA<^^» and DNA<^-'^«^^-«= are consistent with the vahies 
obtained using the restricted fits (typicaUy within 10%). However, the occupancy 
extracted for DNA^-'"^ m die three species sample is smaller (30%), with 
compensating increases in the brightness extracted in R. For the unrestricted 
intermediate-occupancy fits in Table 10, the occupancies for DNA^^-'^ were 
smaller (25%-50%), with compensating increases in the brightness iaR. The 
occupancies extracted using die unrestricted fit from DNA^^-"^ and 
Myisfs. m the fluee-fiiBgment mixture are different from those obtained 

wifli the restricted fit (a 33% decrease and a 160% mcrease). With die unrestricted 
fits, the diffusion times extracted for dnA*^*'"'™^^'^^^ and DNA*^^"^ are similar to 
the theoretical one (500-700 jis vs. 710 ps). The diffusion times for DNA^'^ are 
shorter (-400 jis). The consistency found for brightness is excellent for all fits 
(typically -10%). 

[0141] The effects of triplet state fluctuations and fluorescence saturation [15], 
photobleaching [16, 70], and photo-induced isomerization [67] are weU- 
characterized in FCS; many of these features are evident in the data analyzed here. 
To compare with results using PAID, the data from Tables 8 and 10 were analyzed 
using FCS over the same range of time intervals. modeUng only diffusion (not 
shown). Variations m dif&sion time that coirelated with fliose in Table 8 and. 10 
were found, nidicating that the photophysical properties that affect FCS have similar 
effects on PAID, hnprovement of accuracy and consistency of fitted vahies is 
expected when tiiese effects are mcorporated mto die PAID model. 
[0142] Analysis of RNA polymerase - DNA interactions: PAID was used to 
study the formation of RNA polymerase complexes with DNA, and to show that 
PAID can resolve and quantitate species widi differences in brightness and/or 
diffusion time. DNA was labeled with Cy5 (as the "red" fluorophore r) to yield 
DNA^y^'^B. j^j^^j^^ ^.^ "yellow" fluorophore on the o 

subunit to yield Ra™^. Li the Ra™«-DNA^^' complex, die distance between TMR 
and Cy5 is »100 A, excluding interactions between the fluorophores [57]. We 
examined free DNA*^^ free Ra™». and an interaction mixture fliat contains the 
Ra™^-DNA'^''"» complex along with free species ("RNAP+DNA" sample; Table 
1 1). PAID analysis of DNA'^'^'^b ^covered brightness and difiusion times sunilar 
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to ones of the previous DNA fragments (8.9 kHz vs. 9.4-10.2 kHz, and 440 fis vs. 
^00 ^is, respectively). Analysis of Ra™^ recovered brightness of -10 kHz, while 
the dif&ision time of -710 was larger than that of DNA^^, but smaller than the 
theoretical value (~1 ms; Materials and Methods). In addition to possible causes for 
deviation mentioned in the previous section, thie smaller value may also be due to 
imperfect modeling of the hydrodynamic radius of the protein. 
[0143] The RNAP+DNA sample shows the presence of 4 major species that are also 
observed after resolution of RNAP+DNA samples on polyaciylamide gels (Fig. 12). 
For the KNAP-DNA interaction, distinct species can be identijQed on the basis of 
their gel mobility and their Y:R emission ratio (note that electrophoretic mobility 
depends on size and charge, and is not simply related to free translational diffusion 
measured by PAID). The major species are characterized as i) free DNA; ii) a- 
associated species; iii) KNAP-DNA complex; and iv) aggregates (a heterogeneous 
species with large variety of sizes ranging from smaller oligomers to large 
aggregates). When a 3-component PAID fit was performed on RNAP+DNA sample 
(not shown), a single species with only 7-channel emission was extracted, with 
brightness corresponding to ~2 copies of TMR per species. This result was 
inconsistent with the biochemical system (Fig. 12), with visual inspection of the 
PAID histograms, and with fluorescence-intensity time traces showing few 
exceptionally bright bursts in the 7 channel (such bursts were absent in the buffer or 
DNA^^'^^®, and were extremely rare in the Ra™^ sample). Therefore, a 4-. 
component PAID fit was performed on the RNAP+DNA data, recovering 4 species 
with difiusion times and Y:R brightness ratios that correspond well to gel-resolved 
species of RNAP+DNA samples (Fig. 12). This fit yielded additional infomiation 
inaccessible to ensemble fluorescence imaging or other FFS methods. 
[0144] Species 1 has properties consistent with those of free DNA^^^*^^^ (emission 
only at the R channel, corresponding to a single copy of Cy5); we identify this 
species as "free-DNA". Species 2 shows emission only at the 7 channel 
(corresponding to a single copy of TMR), and difftision time smaller than the 
RNAP-DNA complex or holoenzyme, but similar to DNA (compared with the 600- 
700 us found with Cy3 in previous sections); we identify this species as a 
combination of a "o-associated species" generated during the preparation of RNAP- 
DNA complex (Materials and methods), and of RNAP-DNA complex with non- 
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fluorescent Cy5 (analogous to dNA^^'^'^^'^b ^ p^vious section). Species 
3 shows comparable emission at both 7 and channels (corresponding to a single 
copy of TMR and a single copy of Cy5), and diffiision time (790 ps) slightly larger 
than the free holoenzyme (710 \is) and larger than free DNA (440-480 us with Cy5 
attached, 550-700 \is with Cy3); we identify this species as "RNAP-DNA complex". 
Species 4 shows emission only at the 7 channel, with high brightness (corresponding 
to a mean of ~3 copies of TMR per species), and has the longest average diffusion 
time (850 ^s) of all species; we identify this species as "aggregates". Species 4 is a 
heterogeneous species, likely containing aggregates with varying numbers of 
subunits; different sections of data contained bursts of highly variable height and 
duration (most of the bright bursts from species 4 are of -1 ms duration; however, 
we have observed "yeUow"-labeled species that generated bursts longer that 1 s with 
a peak height of 60 kHz). This variability is manifested in the PAID fits by the 
larger uncertainties found for the brightness and diffusion time of species 4. PAID 
also recovered the occupancy of the various species; however, since PAID 
measurements are performed in a different concentration regime, matrix, and buffer 
than the gel mobility assays, quantitative agreement of PAID results with the gel 
results is not expected. Therefore, this comparison was not performed. 
[0145] The qualitative agreement between PAID and gel-imaging in terais of 
molecular size (measured as PAID-based diffiision time or gel-based electrophoretic 
mobihty) and relative brightness (measured as ratio of PAID-based molecular 
brightness or gel-based fluorescence intensity ratio) of the major species 
demonstrates tiie use of PAID for analysis of complex mixtures. Such a complicated 
mixture will present a great challenge for other popular FFS methods, such as PCS, 
FCCS, PCH/FIDA, 2D-FIDA or FIMDA. These methods are not able to 
simultaneously monitor the diffiision time and brightness on both channels, 
preventing the reliable detection and characterization of the species obtained using 
PAID (2D-FIDA would be able to separate and identify the species, but would lack 
the information on diffiision time). PAID reliably separates the species by 
extracting the relative brightness in the two channels; the diffiision time assists in the 
separation, but also proves very useful in species identification. For example, the 
diffiision time provides additional evidence for the discrimination between the 
aggregates and the dimmer, o-associated species that emit in the same channel. 
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(0146J Fluorescence fluctuation mefliods are valuable tools for in vitro and in vivo 
analysis of macromolecular interactions. The most significant advantage of lliese 
methods is the ability to detect distinct subpopulations and extract their individual 
properties. The subpopulations are best identified and characterized by analysis 
mediods that measure several of each species' properties. PAID is a multi- 
dimensional method that is applicable to muWple channels. The above examples 
demonstrate that PAID provides robust and simultaneous extraction of occupancy, 
diffusion time, and brightness in single- and dual-channel formats for multiple 
species. The examples show that PAID matches or exceeds the statistical accuracy 
of existing methods. 

[0147] In single-channel applications. PAID monitors occupancy, diffusion time, 
and brightness of several species as functions of experimental conditions or time. 
Monomeric species can be distinguished fiom oligomeric species by differences in 
their brightness and difiusion time (Fig. 1 A). The examples on DNA model systems 
demonstrate that brightness measured by fluorescence fluctuation experimaits can 
be used for accurate determination of the number of fluorophores attached to a 
diffusing species, demonstrating the ability to use brightness in fluorescence ■ 
fluctuation methods for extracting stoichiometry. Previously, other groups have 
performed studies that suggest that this is possible, but since no systematic attempts 
to control stoichiometry or purity were undertaken, their results can not be 
considered as rigorous tests of the abiKty of FFS to ectract stoichiometry [24. 25]. 
In contrast the careful design, construction and purification of the DNA model 
systems used in the present examples, combined with PAID analysis establishes this 
capability. The examples also demonstrate that PADD is capable of detecting 
heterogeneity in brightness, while extracting diffusion times. 
[0148] Due to background, due to oligomerization-induced species carrying variable 
numbers of monomers (for example, aggregation induced heterogeneity), and due to 
other comphcations (incomplete labeling, quenching of fluorescence, and 
photobleaching), non-ideal systems contain several species with different brightness. 
To better characterize such systems, the analysis method should detect heterogeneity 
in brightness and diffusion time. Correlations be^veen increases in brightness and 
increases in diffusion time can help verify molecular interactions, providing more 
evidence than if brightness and difftision time were determined individually. FCS is 
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not able to detect heterogeneity in brightness, and relies on using difRision time 
differences to detect subpopulations. PCH/FIDA detect brightness heterogeneity, 
but do not account for difRision (which can differentiate aggregates with quenched 
fluorescence from monomers without quenched fluorescence). The time bin size 
must be kept shorter than the diffusion time, preventing full use of the photon data 
stream. FIMDA has capabilities siroilar to single-channel PAID. PAID and FIMDA 
show similar accuracy for most of the simulations, but FIMDA shows lower 
accuracy for the parameters extracted from simulations in low-occupancy conditions 
(an experimental comparison has not been performed). On the other hand, FIMDA 
aheady accounts for triplet-state fluctuations, whereas these remain to be 
incorporated into the PAID model. It is possible to use a simultaneous fit of FCS 
and FTDA/PCH to extract information about brightness and diffusion time 
simultaneously [31]. However, our fits show that, even if FIDA/PCH fixes the 
occupancy and brightness of two species to their correct values, FCS is not able to 
extract the values for the diffusion times as reliably as PAID or FIMDA (Table 2). 
Using its capabilities, single-channel PAID can be applied to areas such as assembly 
and stoichiometry of membrane proteins (receptors, ion chaimels) [71, 72], and 
amyloid plaque formation [73-75]. FCS has been appHed to the characterization of 
amyloid p-peptide polymerization [76, 77]; a fixed concentration (10 nM) of labeled 
monomer was combined with a varying concentration of unlabeled monomer. The 
presence of aggregates was detected as large changes in diffusion times. Single- 
channel PAID allows both brightness and diffusion time to be quantified, providing 
two observables for the degree of aggregation. 

10149] Dual-channel methods provide improved sensitivity over single-channel ' 
methods for the analysis of the interactions using fluorescence fluctuation methods 
[8]. The examples demonstrate PAID-based extraction of coincidence, diffusion 
time, brightness, and occupancy of several species in a single data set. Control 
experiments were used to restrict many parameters, increasing confidence in the 
remaining fitted parameters. The PAID analysis of the RNAP-DNA interaction 
demonstrates the abihty of PAID to analyze complex systems, detecting 
concentration, brightness and diffusion time for multiple species present in a 
mixture. This ability allows quantitation of the various free and bound species 
present in equifibrium binding reactions, thus paving the way for generating binding 
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constants between interacting partners in a high-throughput, low-volume assay 
format. 

[0150] By extending FCCS to PAID, brightness infonnation can be used to 
discriminate between species and background. In FCCS, leakage and background 
can hamper extraction of the occupancy of complexes [33]; although control 
experiments can extract leakage and background, it is often essential to extract all 
parameters from one data set. requiring control experiments only for checking 
consistency. PAID simultaneously extracts parameters, and is useful in cases where 
precise control e;q>eriments are impractical or impossible. For example, ceUular 
autofluorescence changes as a function of position, preventing background 
measurements for a particular spatial position in a cell that contains different 
fluorescent species. PAID can replace FCCS when it is necessary to extract the 
brightness of a species in addition to occiq)ancy and difiusion time. FCCS was used 
to monitor the endocytosis of Cy2- and Cy5-labeled chblera toxin (CTX) into cells 
[78]. PAID can improve the characterization of flie species by detennining the 
brightness in each channel and the number of CTX subunits per vesicle, while 
analyzing the rare events with long diffusion times and high brightness that baused 
difSculties in the cross-correlation experiments. 

[0151J For systems with more than two channels, PAID can be used by picking a 
series of histograms with channel assignments SIM. Previously, single-pair FRET 
measurements were used to monitor the folded and unfolded subpopulations of the 
protein Caiymotiypsin Inhibitor 2 at differmt concentrations of deaiatmant [69]. By 
splitting the emitted signal by polarisation as weU as by emission spectrum (into 
donor and acceptor; 4 channels total), fhiorescence anisotiopy can be used in 
conjunction with FRET to obtain more accurate distance measurements [36]. Hie 
single-molecule burst identification methods used in these stadies can suffer fiom 
biases due to unequal detection efficiencies for different channels since the burst 
searching algorithms preferentially exclude dimmer species. In extracting 
occupancies, brightness and difiusion time in all four channels, PAID can provide 
unbiased estimates critical for characterization of protein folding. Moreover, while 
single-molecule burst-analysis methods are restricted to low-occupancy samples, 
PAID is applicable to low- and intermediate-occupancy samples, expanding the 
dynamic range of the analysis of interactions. In addition, with appropriate 
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extensions to the PAID model, confonnational dynamics of individual 
subpopulations on time scales shorter than diffusion may be monitored. 
[0152] PAID provides a convenient visual representation of important features of 
diffusing species over a large dynamic range. This is advantageous, since clear 
visual features in histograms provide the first clues to interesting findings (such as 
additional subpopulations) or experimental problems. The PAID histogram focuses 
on photon-rich time intervals, retaining the intuitive nature of correlation functions, 
while simultaneously providing information avaHable fixim photon counting 
histograms. In contrast, the photon counting histograms used in FIDA/PCH and 
FIMDA weigh all time intervals equally; at low occupancy, most of tbe bins 
correspond to time intervals when no molecule is present and few photons are 
counted. In addition, the log axes used in PAID allow for easy visualization of 
effects over a large dynamic range. 

[0153] The model used for the detection volume has a large impact on the accuracy 
of the method. Using our theoretical model for the detection volume rather than a 
Gaussian detection volume improved the considerably. Similar effects ar^ found 
in FCS and other FFS methods [68]. Improvement in the quality of PAID fits will 
come from improved characterization of the detection volume. The FIDA family of 
methods model the detection volume as a volume density per intensity level; a 
polynomial expression is used whose parameters are determined using experiments 
on a standard sample. While attractive for its simplicity, this prevents the modehng 
of possible trajectories through die detection volume. The model used for PAID is 
able to account for different possible trajectories through the detection volume. 
[0154] PAID can impact measurements where simultaneous determination of 
coincidence and diffusion (or other temporal dynamics) are critical. PAID can be 
used to monitor protein-protein interactions, such as protein oligomerization. The 
PAID model can also be extended to include photophysical properties of die dyes 
(triplet-state induced blinking, singlet and triplet state saturation, photobleaching), 
and to incorporate an experimentally-measured detection volume; these 
improvements will increase confidence and might uncover additional dynamics or 
species not assumed in the fit. The model can be extended to account for two 
monitor channels, and to become compatible with studies of immobilized molecules 
and systems mvolving flow. These extensions of PAID will aUow in vitro analysis 
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of dynamics (such as aspects of protein folding), protein-protein interactions and 
protein-DNA interactions, as weU as analysis of interactions in a ceUuIar 
environnient. 

[0155] "Transit time" for tbe purposes of this specification includes difiusion time 
and/or flow time of species within the detection volume. In addition, the term 
"flnorophore" is intended to cover not only organic fluorescent fluorophores. but 
also any label or tag, intrinsic or extrinsic, that emits or scatters photons when 
subjected to light excitation. The PAID method is also useful for measuring 
temporal dynamics of species within the detection volume. ' ' • 
[0156] Having thus described exenq)laiy embodiments of the present invention, it 
should be noted by those skilled in the art that the within disclosures are exemplary 
only and that various other alternatives, adaptations and modifications may be made 
within the scope of the present invention. Accordingly, the present invention is not 
limited to the above preferred embodiments and exanqiles, but only limited by fhe 
following cleiims. 
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